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I  Introduction 


Reported  below  are  the  results  of  a  study  carried  out  at  the  Center  for  Communication  and  Signal  Processing 
Research  at  New  Jersey  Institute  of  Technology  between  July  1994  and  December  1995.  This  research  can  be 
considered  a  continuation  of  previous  research  performed  during  1991/93  [1]  ,  1992/1993  [2]  and  1993/1994 
[3].  Its  aim  is  to  further  widen  the  investigation  of  methods  for  separation  of  multiuser  signals.  This  can 
occur  in  satellite  communication  microwave  terrestrial  links  or  in  wireless  multiuser  channels.  Although 
some  emphasis  was  placed  on  using  code  division  multiple  access  (CDMA),  this  is  not  a  requirement  in  the 

proposed  approaches. 

Historically  the  group  at  the  Center  has  been  involved  in  investigating  algorithms  for  interference  cancela¬ 
tion  and  signal  separation.  The  algorithm  termed  “Bootstrap  Algorithm”  was  first  proposed  by  the  principle 
investigator  in  1981  [4],  and  later  it  was  used  for  cancelling  cross  polarization  in  satellite  communications  [5] 
and  in  microwave  terrestrial  radio  link  [6,  7,  8], 

During  the  91/92  and  92/93  phases  of  this  research,  a  thorough  study  of  the  “bootstrap  algorithm”  was 
performed.  Three  different  structures  were  proposed  in  [1];  namely,  the  Backward/Backward  (B/B),  the 
Forward/Forward  (F/F)  and  the  Forward/Backward  (F/B)  structure.  They  are  depicted  in  Figures  1,  2  and 
3,  as  they  applied  to  “dually  polarized  signals.”  The  bootstrap  separator  of  co-channel  signals  is  principly 
composed  of  cancellers,  each  one  using  the  outputs  of  the  other  cancellers  to  facilitate  control  of  its  adaptive 
weights.  In  fact  such  structures  performs  as  a  “signal  separator”  rather  than  an  interference  canceler.  Since 
for  its  operation,  there  is  no  need  for  a  reference  signal  (as  in  the  LMS  canceler),  it  is  sometimes  justifiably 
called  a  blind  separator.  Note  that  in  the  B/B  algorithm,  the  weights  are  controlled  by  power  minimization; 
in  the  F/F  algorithm  they  are  controlled  by  decorrelating  the  outputs;  while  in  F/B  by  a  combination  of 
both.  It  was  shown  that  for  these  algorithms  to  converge  to  a  state  of  signal  separation  a  signal  distinguisher, 
termed  discriminator  (DIS),  is  needed.  Such  a  discriminator  uses  some  known,  simple  estimates  of  the  signals 
themselves  as  in  the  case  of  the  regular  LMS  noise  canceler.  In  general  the  adaptive  weights  are  controlled 
with  recursive  equations  as  follows; 
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for  the  case  of  the  power  criterion 


j^i  z,  j  =  1 . . .  K, 


and  for  the  decorrelator  control 


%3 


+  f^f 


\iDzi)zji 


j  i 


(1) 


where  Zi  i  —  I  -  ■  -  K  {K  is  the  number  of  co-channels)  is  the  output  of  the  decorrelator,  and  D  is  the 
discriminator  operator.  In  the  work  reported  in  [1]  and  [2]  the  discrimihator  was  assumed  heuristically  as 
imposing  a  slightly  different  gain  on  the  signal  and  interference  at  each  output  (also  see  published  references 
[9,  10,  11,  12]).  There  it  was  possible  to  examine  convergence  behavior  and  error  probability  performance 
and  compare  it  to  other  separators  known  in  the  literature.  In  [13]  a  third-order  function  was  used  as  a 
discriminator.  Some  preliminary  results  were  also  obtained,  which  describe  the  92/93  phase  of  the  research, 
in  applying  these  algorithms  for  separating  co-channel  CDMA  signals,  and  blind  channel  equalization. 

During  93/94  [3]  further  research  was  performed  in  two  directions  of  application:  separation  of  a  multiuser 
CDMA  signal  and  blind  equalization  of  dispersive  communication  channels.  For  the  first  topic,  we  proposed 
a  two  stage  canceler  for  synchronous  multiuser  CDMA  [14].  The  first  stage  is  a  decorrelating  detector  while 
the  second  is  an  adaptive  interference  canceler.  To  further  improve  the  performance,  a  soft  decision  was 
used  instead  of  the  hard  decision  that  precedes  the  second  stage  [15].  The  slope  of  the  soft  decision  was 
chosen  heuristically  based  on  measurements  performanced  on  the  output  of  the  first  stage.  The  bootstrap 
separator  was  also  used  as  a  second  stage  instead  of  a  minimum  power  canceler.  The  work  also  extended  to 
the  asynchronous  case.  A  one  shot  matched  filter  was  used  to  model  this  case. 

For  blind  equalization,  we  proposed  a  blind  decision  feedback  equalizer  base  on  the  decorrelator  [16], 
and  we  introduced  a  blind  maximum  likelihood  sequence  estimator  [l'7],'and  a  form  of  reduced  complexity 
sequence  estimator  that  uses  state  partitioning  [18]. 

During  this  phase  of  research  we  concentrate  on  three  main  topics; 


1.  We  perform  a  thorough  study  of  the  F/F  bootstrap  which  uses  the  signum  function  as  a  discrimina¬ 
tor.  We  examine  its  error  performance  when  it  is  applied  to  multiuser  CDMA,  and  we  compare  its 


performance  to  what  is  known  as  a  correlating  detector  on  one  the  hand  and  to  the  MMSE  detector 
on  the  other.  We  show  that  in  contrast  to  the  other  detector,  the  bootstrap  can  perform  better  when 
the  channel  cross-correlation  matrix  is  not  accurately  known.  We  further  suggest  a  method  for  bit 
combining  estimated  data  in  an  asynchronous  channel,  which  results  in  an  improvement  m  output 
error  probability  of  the  decorrelating  detector.  A  method  for  signal  separation  that  is  closely  related 
to  the  F/F  bootstrap  stage  algorithm,  termed  gradient-based  decorrelating  algorithm,  was  also  studied 
and  used  for  multiuser  CDMA.  It  is  shown  that  such  an  approach  depicts  a  smaller  eigenvalue-spread 
than  that  of  the  corresponding  method. 

2.  An  adaptive  array  is  also  used  for  performance  improvement  in  multiuser  communication. 

3.  A  new  approach  is  suggested  to  handle  multiuser  signal  separation  in  asynchronous  CDMA.  It  can,  in 
particular,  handle  singular  cases. 
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II  Technical  Summary 

11*1  The  FF  Bootstrap  Algorithm  With  Signum  Function  Discriminator 

For  the  general  multi-user  the  equivalent  low-pass  signal,  the  input  of  the  matched  filter,  is  given  by 

K 

(0  =  (0  {i)sk  {t  +  iT-  n)  +n(t),  (2) 

fc  =  l  i 

where  K  is  the  number  of  users,  and  a*,  6*,,  sj  and  r*  are  the  signal  amplitude,  user  bit,  signature  waveform, 
and  relative  delay  of  the  user,  n  (i)  is  the  zero  mean  AWGN,  with  a  two-sided  power  spectral  density  of 
No/2. 

For  the  synchronous  channel  encountered  in  down-link  communication  channels  the  r*  are  zero  for  all  k. 
In  the  asynchronous  up-link  channels  the  Tk  are  not  necessarily  equal.  For  the  former  case  the  output  of  the 
ir**  matched  filter  is  the  composite  of  bit  bk  and  all  interfering  bits  given  by  a  linear  combination  through 
the  cross-correlation  factor  pkj  •  In  matrix  notation  we  can  write 

x(j)  =  PA6(i) -f  n(j')  ,  (3) 


1 

Pl2 

Pl3 

PlK 

Pl2 

1 

P23 

P2K 

Pl3 

P23 

PK-1,K 

.  PlK 

P2K 

Pk-i,k 

1 

and  A  =  diag  (^37 ,  •  •  ■ ,  ,  6  (i)  =  [bi  (i)  ,  •  ■  • ,  b^  (i)]  and  n  (i)  is  a  zero-mean  Gaussian  noise  vector  with 

covariance  Rn  =  Pa^.  For  the  asynchronous  case  one  may  use  the  one-shot  matched  filter  proposed  by  Verdii 
in  [19]  .  In  this  case  (3)  will  contain  a  linear  combination  of  bk  (i),  bj  (i)  and  bj  (i  -  1)  for  y  =  1 . . .  K,  j  ^  k. 
That  is,  the  matrix  P  will  be  replaced  by  another  square  matrix  with  l-f2  (if  —  1)  =  2A'  —  1  dimensions.  The 
output  of  this  one  shot  matched  filter  Xk  (i)  will  contain  a  linear  combination  through  this  (2if  -  1)  x  (2A'  -  1) 
cross-correlation  matrix  of  bk  (*)  and  bits  i  and  i  -  1  of  all  other  users.  If  the  r,  are  known  then  all  the 
elements  of  this  matrix  are  known  provided  one  knows  the  codes.  Therefore,  for  both  synchronous  and 

asynchronous  cases  the  output  of  the  matched  filters  is  given  by  (3),  with  the  corresponding  definition  of 

b(i). 
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II.l.l  The  principle  bootstrap  structure 

Following  the  matched  filter  we  have  (see  Fig.  4) 

x  =  Vx  =  V  PAb+Vn=VPAb  +  <^.  (5) 

We  use  y  =  7  —  VF  with 


= 


0  Wi2  U)13 

ti)21  0  ^^)23 


WlK 

W2K 


(6) 


:  0  wk-\,k 

wki  wk,k-i  0 

For  the  two-user  case  the  principle  structure  using  the  signum  discriminator  in  given  in  Fig.  5.  That  is,  the 

weights  are  controlled  by  the  recursion 


«^l2(n  +  l)  =  i«l2(”)  -  P^isgn(z2) 

tt2i(n  +  l)  =  iii2i(n)  - /iZ2Sgn(/i). 


(7) 


The  mean  steady  state  value  of  these  weights  are  given  by 


E(^isgn(ij))  =  0  i  i  =  1, 2.  (8) 

In  Appendix  A  [20]  it  is  shown  that  when  zj  has  a  high  SNIR  then  sgn(zy)  ~  hj  and  the  corresponding 
steady  state  weights  equal  the  corresponding  cross-coupling  (wi2  =  tD2i  =  p),  in  which  case  total  cancellation 
of  interference  occurs.  This  is  sometimes  called  near-far  resistance,  since  the  probability  of  error  does  not 
depend  on  the  level  of  interference.  For  the  general  multiuser  case  we  have 

SNRk-P^^Pk,  (9) 

where  the  index  k  for  a  matrix  or  vector  means  the  ^-th  row  and  column  k  are  taken  out  of  the  matrix,  or 
that  the  k-th  element  is  taken  out  of  the  vector. 

At  this  limiting  condition  or  high  SNIR  we  have  for  the  signal-to-noise  ratio 

SNR,  =  ‘A'rAPpR,) 


9 


On  the  other  hand  when  Zj  does  not  have  a  large  SNIR  to  justify  sgn(zj)  =  hj  then  (8)  is  not  valid  in  steady 


For  the  two-user  case,  let  sgn(iri)  =  hi,  however  sgn(2:2)  62  then  E(0isgn(z2))  =  0  implies  W21  =  p.  But 

E(^2Sgn(zi))  =  0  will  only  imply  u;i2  =  p  -f-  6  (u;i2  p)-  Performing  the  required  analysis  in  Appendix  A, 


0  =  — - ; - : - _ 

1  +  (1  -  2Q  (yZSNR^'j  e' 

where  LSNR  is  the  limiting-signal-to-noise  ratio  at  the  output: 


LSNR 


02  (1  -  P^) 


From  (11)  we  notice  that  6  has  the  opposite  sign  of  p  and  that  its  absolute  value  |^|  decreases  from  \p\ 
to  zero  when  the  LSNR  increases  from  zero  to  infinity. 

The  corresponding  probability  of  error 

-2^ 

Ip  (  i'i- -  -  p6)  -  y/^6\ 

2  y  <r^^/i  —  p2  q.  ^2  j  (  ) 

For  small  02,  6  —*  —p.  Then  Pei  — Q  the  error  corresponding  to  a  simple  BPSK  signal.  For 

large  02,  (5-^0,  and  Pg^  ^  Q  -  p^/crj,  the  error  corresponding  to  the  limiting  case  (Verdii’s 

decorrelating  detector),  see  Fig.  6. 

Note  that  if  one  knows  the  matrix  P,  then  instead  of  V,  we  may  use  p-^  (this  is  the  Verdudecorrelating 
detector)  whose  probability  of  error  is  constant  and  equals  the  limiting  value.  Hence  we  conclude  that  even 
if  P  IS  known  the  bootstrap  separator  outperforms  Verdu’s  detector  for  almost  all  SNR.  A  thorough  study 
of  this  bootstrap  structure,  including  the  multiuser  case,  is  given  in  Appendix  A. 


II.1.2  Comparative  Study  of  Linear  Minimum  Mean  Square  Error  (LMMSE)  and  the  F/F 
Bootstrap  Algorithm  with  Signum  Discriminator 

The  LMMSE  multiuser  detector  decodes  h  based  on 


b  —  sgn  (Olmmse^  , 
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where 


Olmmse  —  Wx  —  {P  +  cr^ A  x. 


It’s  probability  of  error  Pe{k)  has  been  derived  [21],  (see  also  appendix  B); 

p-m  =  2^  S  ■  (15) 

Under  weak  and  strong  interference  conditions,  the  above  error  probability  expression  will  reduce  into  the 
single-user  BPSK  limit  or  into  that  of  the  decorrelating  detector,  respectively. 

f  Q  (■^^)  >  weak  interference 


Pe{k)  = 


—  P k  Ph 


strong  interference 


In  Fig.  7  the  performance  of  the  LMMSE  detector  as  compared  to  the  bootstrap  detector  is  shown. 


II.1.3  The  Bootstrap  Algorithm  for  an  Unknown  Channel 

In  many  applications  (not  necessarly  CDMA)  the  matrix  P  is  unknown.  In  CDMA  applications  the  matrix 
P  might  not  be  accurately  known,  as  a  result  of  the  unknown  effect  of  the  fading  channel  or  as  a  result 
of  errors  in  estimating  the  relative  delays  between  users  in  the  asynchronous  uplink  mobile  communication 
channel. 

To  see  the  advantage  of  the  bootstrap  algorithm  in  these  situations  we  first  write  (2)  for  the  two-user 
case  as 
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As  shown  in  Fig.  8,  we  apply  r{t)  to  a  bank  of  filters  by  using  as  a  second  input  S2' (t),  S2' (t), 
respectively,  the  signature  function  of  Eq.  (18)  extended  to  instead  of  T2,  normalized  by  instead  of  £2- 
As  a  result  we  have  at  the  output  of  the  matched  filters  for  the  zero-th  bit 


x'(0)  =  P^Ab(0)  +  n(0), 


(19) 


where 


A  =  diag[i/ar,  s/a^,  ^02(1  -  £2)]  and 


1  P21 

P21  P22^ 

.  P12  P22^ 


Pl2 

P22^  J 


(20) 


b(0)  =  [61(0),  62(~1)i  ^2(0)]^  .  n(t)  is  zero-mean  Gaussian  vector  with  covariance 


R 


n 


1  P21 

P21  1 

Pi  2  0 


2  ■ 


(21) 


If  the  relative  delay  T2  is  known  then  we  choose  =  r2,  and  as  a  result  p'21  =  P21 ,  P12  =  Pi2>  P22^  =  P22^  =  1> 
P22^  —  P22^  =  0  and  £^  =  £2  =  1.  The  matrix  Pi>  is  known.  If  ]p2i|  +  IP12I  <  1  then  Pi?  is  invertible.  A 


simple  linear  transformation  by  P^^  results  in 


z(0)  =  Ab(0)-1-P^^n(0), 


(22) 


where  Pd  is  the  modified  cross-coupling  matrix.  Since  the  components  of  b(0)  are  assumed  uncorrelated 
then  the  components  of  z(0)  are  also  so.  A  detailed  analysis  of  applying  partial  codes  corresponding  to  r' 
instead  of  r  as  given  in  Fig.  8  is  in  Appendix  C.  (see  also  [22]). 

The  results  of  using  the  bootstrap  when  <  T2  so  that  £3  =  0.35  instead  of  £2  =  .4  is  given  in  Fig.  9 
and  Fig.  10.  Note  that  when  using  the  conventional  detector  with  P'~^  (instead  of  P~^)  the  probability 
of  error  is  approximately  0.5,  while  with  the  bootstrap  bit  error  rate  is  adequate  and  almost  equal  to  that 
when  P  is  known. 
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11.1.4  Maximum  Ratio  Data  Combining  for  Performance  Improvement  of  One-shot  Asyn¬ 
chronous  Detectors 

Using  a  one-shot  matched  filter  followed  by  a  decorrelating  detector  results  in  outputs  corresponding  to  the 
main  user  and  the  left  (earlier  bit)  and  right  (current  bit)  portion  of  the  code  of  the  other  users.  Due  to 
reduced  power,  these  outputs  will  result  in  a  degradated  probability  of  error.  In  Appendix  D  we  propose  a 
method  of  combining  the  data  of  the  split  code  users.  We  suggest  using  a  combining  gain  in  proportion  to 
the  relative  user  delays.  The  resulting  decorrelating  detector  always  performs  better  than  the  split-output 
detector.  For  the  two  user  asynchronous  case,  following  the  matched  filter  we  have 

x(0)  =  PAb(0)-l-n(0),  (23) 


where 


1  0 
0  1 


(24) 


A  =  di*g[»„ <.,,«]  =  and  b(0)  =  lhmM-l).h{0)r ^  The  „(<)  ie  a 


zero-mean  Gaussian  vector  with  covariance,  PiVo/2. 

T 

=  4=  / 

V^2  Jo 

T 

pf  =  -7f=  /  (23) 

Vl  -  ^2  Jo 

If  we  use  the  bootstrap  with  V  =  J-W  then  following  the  result  of  eq.  (9)  we  have 

zi(0)  =  ^{l-pJPT"pi)bii0)  +  nDi{0) 

Z2(0)  =  V2)^2(-1)  +  nD2(0)  (26) 


Z3(0)  =  \/a2(l  -  e2)(l  -  pJPs  V3)^2(0)  +  nc3(0), 


with  the  corresponding  SNR 


SNRfc 


No/2 


k  =  1,2,3. 


We  suggest  combining  the  data  at  the  output  of  the  bootstrap  as  in  Fig.  11 


22c(0)  —  2:2(0)  +  72:3(— !)• 


(27) 


(28) 
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In  Appendix  D  (see  also  [23])  it  is  shown  that  choosing  7  =  will  result  in  maximum  SNRzc  It  is  also 

shown  that  such  a  choice  will  make  SNR2C  greater  than  the  uncombined  SNR2  or  SNR3  for  any  e  and  any 
P .  As  depicted  in  Fig.  12,  combining  data  gives  better  error  probability  performance. 

11.1.5  Adaptive  Setting  of  Threshold  for  Soft  Tentative  Decision  of  a  Multistage  Canceler 

In  [3]  we  presented  a  multistage  separator  of  a  correlation  detector  (or  bootstrap  separator)  followed  by 
a  hard  limiter  tentative  decision  followed  by  adaptive  canceler.  We  also  show  that  using  a  soft  limiter 
instead,  the  performance  would  improve.  The  soft  limiter  used  some  heuristic  threshold.  In  this  phase 
of  research  (see  Appendix  E),  we  propse  to  adjust  the  threshold  level  so  as  to  optimize  the  performance. 
The  structure  of  two-stage  (decorrelator/canceler)  with  adaptive  soft-limiter  is  depicted  in  Fig.  13  with  the 
adjustable  threshold  shown  in  Fig.  14.  The  optimal  canceler  weights  and  optimum  threshold  values  were 
found  analytically.  In  simulation  the  following  two  steepest  descent  algorithms  were  used. 

=  <(Af,.)  +  //V..,F^{2/?(Af,.,u;,^.)} 

+  n  Va,,-  (29) 

Comparing  the  performance  when  using  a  hard  limiter,  the  heuristic  threshold  and  the  optimal  threshold 
are  given  in  Fig.  15. 

11. 1.6  Convergence  of  a  Stochastic  Gradient-Based  Decorrelation  Algorithm 

We  analyzed  the  performance  of  a  stochastic  gradient-based  decorrelation  algorithm  for  blind  signal  separa¬ 
tion  [25,  24]  (see  also  Appendix  F).  The  algorithm  has  a  computational  complexity  comparable  with  the  LMS 
algorithm,  but  is  shown  to  be  faster  than  the  latter.  The  different  speeds  of  convergence  of  the  decorrelation 
algorithm  and  the  LMS  are  shown  to  be  a  consequence  of  the  different  eigenvalue  spreads  associated  with 
the  implementation  of  the  two  algorithms.  While  the  problem  statement  and  proofs  are  quite  general,  the 
specific  application  considered  is  a  synchronous  multiuser  multicarrier  CDMA  communication  system  in  a 
fading  environment.  In  this  application  the  decorrelator  is  used  as  a  multiuser  detector.  The  convergence 
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speed  of  the  two  algorithms  is  illustrated  in  Fig.  16,  which  shows  the  average  learning  curves  of  the  proba¬ 
bility  of  error  for  SNR  =  8  dB  and  SSR  =  -5  dB,  where  SSR  is  the  ratio  of  the  bit  energy  of  the  designated 
user  to  the  user  who  serves  as  interference.  The  higher  rate  of  convergence  of  the  decorrelator  is  evident. 

II.2  Adaptive  Antenna  for  Performance  Inprovement  of  Multiuser  Communi¬ 
cations 


Antenna  arrays  provide  diversity  paths  to  combat  multipath  fading  and  are  capable  of  suppressing  spatial 
interference  sources.  Under  consideration  is  a  narrowband  multiple  access  system  such  as  TDM  A.  The 
channel  is  assumed  to  be  flat  fading.  In  [26]  (see  also  Appendix  G)  we  proposed  an  eigenanalysis-based 
filter,  referred  to  as  an  eigencanceler,  as  an  alternative  to  the  Sample  Matrix  Inversion  (SMI)  method  for 
small  data  sets.  In  the  paper  we  develop  analytical  expressions  for  the  average  probability  of  error  of  the  SMI 
and  the  eigencanceler  methods  for  the  case  of  single  interference  and  a  flat  slowly  fading  Rayleigh  channel. 
Consider  the  (I-l-N)  correlation  matrix  Rn  When  the  number  of  interferences  L  is  less  than  A,  the  number 
of  antennas  in  the  array,  the  propagation  vectors  uj ,  j  =  1, ...,  T,  are  linearly  independent  (with  probability 
one)  and  span  a  signal  subspace,  referred  to  as  the  interference  subspace.  Let  Xi  >...>  Xl  >  Xl+i  = 
...  =  Ajv  =  0-2  be  the  eigenvalues  of  Rn  The  L  largest  Ai  are  referred  to  as  principal  eigenvalues  and 
their  associated  eigenvectors  are  referred  to  as  principal  eigenvectors.  The  interference  subspace  is  spanned 
by  the  principal  eigenvectors.  The  noise  subspace  is  defined  as  the  signal  subspace  associated  with  equal 
(jV’_£,q-l)  eigenvalues.  Let  Q^.  and  be  matrices  whose  columns  are  the  dominant  and  noise  eigenvectors, 

respectively.  Then  in  [26]  the  eigencanceler  is  defined  as  the  optimal  weight  vector  constrained  to  the  noise 
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The  random  variable  p  is  conditioned  on  both  the  channel  and  the  covariance  matrix  estimate,  and  is  bounded 
0  <  p  <  1.  The  normalized  SNIR  for  the  eigenanalysis-based  method  is  given  by  the  expression 

U? QvQv  ^QvQv 

The  density  of  p  for  the  eigenanalysis-based  method  is  given  by 

/p(p|7)  =  ^(1  -  (33) 

where  r(.)  is  the  standard  gamma  function.  Since  -y  =  py  and  fy  (7)  and  fp  (p)  are  known,  the  density  of  7 
can  be  found  from  the  expression 

pOO  1  ^ 

fy(y)=  /  -/p(^  It)  A  (7)^7-  (34) 

Jo  1  1 

Combining  eq.  (34)  and  the  conditioned  probability  of  error  Pe(e  I  7)>  the  average  probability  of  error  for 
the  SMI  method  can  be  written  as 

=  r  r  I  7)/p(  J|7)rf7]rf7.  (35) 

Jo  1  Jo  7 

Then,  the  average  probability  of  error  for  the  eigenanalysis-based  method  can  be  found  from  expression  (35), 
with  A(^|7)  =  /p(p|7)  given  by  eq.  (33). 

The  match  between  theory  and  simulations  is  illustrated  by  Figure  17.  The  curves  shown  represent  the 
theoretical  probability  of  error,  and  the  probability  of  error  based  on  simnlations  for  TV  =  9,  RT  =  20,  and 
INR  =  2  dB. 

II.3  Multi-shot  Approach  for  Separating  Asynchronous  Multiuser  CDMA  Sig¬ 
nals 

It  has  been  noticed  that  the  multiuser  detector  proposed  in  [22]  to  handle  asynchronous  transmission  might 
run  into  the  singularity  problem  in  its  enlarged  “correlation”  matrix  for  some  relative  delays  between  multiple 
access  users.  In  order  to  mitigate  the  singularity  problem,  we  proposed  a  new  approach,  termed  multi-shot 
approach,  to  separate  and  detect  multiuser  signals  in  asynchronous  CDMA  Communication  Systems.  The 
multi-shot  approach  first  explores  the  structure  inherent  in  the  matrix  decomposition  of  the  properly  arranged 
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data,  obtain  from  multi-shot  matched  filtering.  Then  we  make  a  matrix  approximation  to  make  the  tentative 
multiuser  information  bit  estimation  tractable  and  simple.  After  analyzing  the  cuase  of  bias  introduced  by 
matrix  approximation,  we  further  modified  the  filtering  matrix  and  reduced  the  bias  effect  and  improved 
the  overall  detection  performance  yet  still  minimized  the  computational  complexity  (linear  to  the  number  of 
users) . 

The  system  diagram  is  shown  in  figure  5  of  Appendix  G  for  the  simpliest  case  which  use  the  minimum 
number  of  consecutive  symbols  (N=3).  Notice  that  all  the  matrix  filterings  involved  in  the  system  are  only 
of  dimenstion  K  (number  of  users).  This  simplification  in  system  implementation  is  achieved  by  exploring 
the  block-tridiagonal  structure  of  the  original  matrix  P  of  a  larger  dimension. 
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Ill  Conclusions  and  Recommendations 


During  this  period  we  made  significant  progress  in  applying  the  bootstrap  algorithm  to  multiuser  CDMA 
communications.  In  particular,  we  use  the  F/F  structure  and  use  the  signum  function  as  a  discriminator. 
The  simplicity  and  superiority  of  the  algorithm  was  emphasized  when  we  compared  its  performance  to  the 
MMSE  on  one  hand  and  to  Verdii’s  decorrelator  on  the  other. 

The  adequacy  of  the  algorithm  to  operate  in  an  unknown  channel  was  also  demonstrated.  Some  ideas 
for  improving  error  probability  in  the  asynchronous  channel  were  also  suggested. 

In  addition,  two  new  topics  not  directly  related  to  the  bootstrap,  but  important  research  initiatives, 
were  studied.  The  first  concerns  with  using  adaptive  arrays  to  improve  separator  performance.  The  second, 
termed  multi-shot  separator,  is  a  novel  idea  that  can  be  implemented  for  an  asynchronous  multiuser  channel, 
and  can  in,  particular,  handle  cases  wherein  the  cross  correlation  matrix  is  singular. 

In  summary,  during  this  period  of  research,  we  obtained  an  important  result  in  applying  the  bootstrap 
algorithm  for  multiuser  CDMA  and  initiated  new  work  on  other  topics  related  to  separation  of  signals. 

A  list  of  publications  resulting  from  these  results  is  given  in  section  V.  For  further  study  we  recommend 
the  following: 

•  Continue  examining  the  application  of  antenna  arrays  to  improve  separator  performance. 

•  Further  examine  the  performance  of  the  multi-shot  separator. 

•  Examine  the  effect  of  channel  fading  and  multipath  on  separator  performance. 

•  Examine  the  effect  of  uncertainty  in  relative  delay  on  performance  of  the  bootstrap  algorithm. 
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Figure  1;  B/B  dual  polarization  separator 


Figure  2:  F/F  dual  polarization  separator 


Figure  3:  F/B  dual  polarization  separator 
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(b).  Medium  capacity  case 


Figure  7:  Performance  of  the  LMMSE  detector  and  the  adaptive  bootstrap 
detector  in  a  multi-user  CDMA  system.  Also  shown  in  the  figures 
are  the  results  of  the  decorrelating  detector  and  the  BPSK  limit 
Parameters  used:  SNRi  =  8dB,  AT  =  3,  M=  1000000  independent 
trials.  26 


SNR2-SNm 

Figure  9:  Performance  Comparison  of  Adaptive  and  Con¬ 
ventional  Detectors  for  fixed  SNRl-variabe  inter¬ 
ference  (Unknown-delay  case) 


Figure  10:  Performance  Comparison  of  Adaptive  and  Con¬ 
ventional  Detectors  for  fixed  SNRl-variabe  inter¬ 
ference  (Unknown-delay  case) 
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Figure  14:  Soft  Limiter  with  controlled  threshold 


Figure  15:  Output  Energy  of  desired  user,  E{2/?}  with 
Hard  Limiter,  Heuristic  Threshold  and  Optimal 
Threshold;  SNRi  =  8  dB,  p  =  0.7 


Figure  16:  Learning  curve  of  the  probability  of  error  of  the 
first  user  for  N=4  users,  SNR=8dB,  SSR=-5dB, 


BIT  ERROR  RATE 


Figure  17;  The  average  BER  versus  the  average  received 
SNR  for  adaptive  array  with  one  interference  with 
K=20,  INR=2dB.  (analytical  results) 
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Abstract 

Decorrelating  detectors  were  proposed  as  tools  to  combat  the  near-far 
problem  in  multi-user  CDMA  communication  systems.  They  were  also 
used  as  the  first  stage  of  many  two-stage,  multi-user  signal  separators  to 
facilitate  tentative  decisions  for  the  canceDer  stage.  These  decorrelators 
result  in  a  complete  separation  of  signals,  but  also  in  relatively  high  ad¬ 
ditive  noise.  An  improved  multi-user  detector  is  proposed  in  this  paper. 
It  may  not  totally  reject  other  signals,  but  it  does  result  in  a  better  er¬ 
ror  probability,  particularly  in  regions  where  interference  energies  are  low. 
The  proposed  improved  decorrelators  can  be  implemented  by  an  adaptive 
algorithm  and  hence  are  better  suited  for  a  changing  environment. 


1  Introduction 

A  conventional  single-user  detector  for  the  Code  Division  Multiple  Access  (CDMA) 
scheme  consists  of  a  bank  of  filters,  each  one  matched  to  the  signature  sequence 
of  a  particular  user.  Besides  the  desired  signal,  the  sampled  output  of  each 
matched  filter  contains  the  residual  interference  from  all  other  users.  The  pres¬ 
ence  of  strong  interference  often  makes  it  impossible  to  detect  weak  users,  a 
condition  referred  to  as  the  near-far  problem.  To  eliminate  this  problem  Verdii 
proposed  a  receiver  [1]  that  is  optimum  in  the  multi-user  environment.  It  is 
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optimum  in  the  sense  that  it  completely  separates  the  signals  from  different 
users’  signals.  Achievement  of  performance  with  this  optimal  receiver  comes  at 
the  expense  of  high  computational  complexity.  A  class  of  sub-optimum  receivers 
that  use  decorrelating  detectors  and  which  are  based  on  linear  transformation  of 
the  sampled  matched  filters’  output  was  considered  in  [2]  and  [3] .  These  decor- 
relating  detectors  were  also  used  as  a  first  stage  of  many  two-stage,  multi-user 
detectors  to  facilitate  tentative  decisions  for  the  second  (canceller)  stage.  This 
second  stage  can  be  a  fixed  canceller,  assuming  knowledge  of  signal  energies  ([4], 
[5]),  an  adaptive  canceller  assuming  no  knowledge  of  energies  ([6]),  an  adaptive 
canceller  with  soft-decision  ([7])  or  a  canceller  with  optimum  soft-decision  ([8]). 
Verdii  s  decorrelating  detector  uses  the  inverse  of  the  cross-correlation  matrix  to 
separate  signals  from  the  different  users.  However,  separation  can  be  achieved 
whenever  the  transformation  diagonalises  the  cross-correlation  matrix.  We  will 
discuss  a  different  transformation  that  will  also  separate  the  signals  and  show 
that  the  output  SNR  and  hence  the  error  probability  remains  unchanged. 

However,  we  are  mainly  interested  in  a  low  error  probability,  not  so  much  in 
total  signal  separation.  Therefore,  the  resulting  signal-to-ii^terference-plus-noise 
ratio  is  the  parameter  that  should  be  minimised.  To  do  this  we  suggest  a  linear 
transformation  that  minimises  the  correlation  between  the  different  transform 
outputs  and  the  hard  limiter  outputs.  This  transformation  is  shown  to  perform 
better  than  the  optimal  separation  decorrelator,  particularly  for  high  signal-to- 
interference  ratios. 


2  Matched  Filters 

For  the  general  multi-user  the  equivalent  low-pass  signal  at  the  input  of  the 
matched  filter  is  given  by 

K 

^  X]  (*)«*  {t  +  iT-Tk)  +  n  (t)  (1) 

i 

where  K  is  the  number  of  users,  and  aj,,  6*,,  sj,  and  r*  are  the  signal  amplitude, 
user  bit,  signature  waveform  and  relative  delay  of  the  user,  n  {t)  is  the  zero 
mean  AWGN,  with  a  two-sided  power  spectral  density  of  Ao/2. 

For  the  synchronous  channel  encountered  in  down-link  communication  chan¬ 
nels  the  Tk  are  zero  for  all  k.  In  the  asynchronous  up-link  channels  the  r*  are 
not  necessarily  equal.  For  the  former  case  the  output  of  the  Ar*'*  matched  filter 
is  the  composite  of  bit  hk  and  all  interfering  bits  given  by  a  linear  combination 
through  the  cross-correlation  factor  pkj.  In  matrix  notation  we  can  write: 

x{i)  =  PAb{i)  +  n{i)  (2) 
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where 

1  Pi2  Pi3  •  •  ■  pm 

Pi2  1  P23  ••  •  Pm 

P  —  Pl3  P23 

;  ;  PK-i,K 

PlK  P2K  PK-1,K  1 

and  A  =  ^>(0  =  [&i  (0  - ' ' '  -  MO]  and  n  (i)  is  a  zero- 

mean  Gaussian  noise  vector  with  covariance  Rn  =  P(t  .  For  the  a.synchronous 
case  one  may  use  the  one-shot  matched  filter  proposed  by  Verdii  in  [9]  .  In 
this  case  (1)  will  contain  a  linear  combination  of  bk  (0,  bj  (i)  and  bj  (z  -  1) 
fQj.  j  k.  That  is,  the  matrix' F  will  be  replaced  by  another 

square  matrix  with  1  -t-  2  (K  -  1)  =  2/^  -  1  dimensions.  The  output  of  this 
one  shot  matched  filter  Xk{i)  will  contain  a  linear  combination  through  this 
(2K  -  1)  X  -  1)  cross  correlation  matrix  of  bk  {i)  and  bits  i  and  z  -  1  of 
all  other  users.  If  the  Tk  are  known  then  all  the  elements  of  this  matrix  are 
known  provided  one  knows  the  codes.  Therefore,  for  both  synchronous  and 
asynchronous  cases  the  output  of  the  matched  filters  is  given  by  (2)  with  the 
corresponding  definition  of  b(i). 


3  Multi-user  Decorrelating  Detectors 


3.1  Verdu’s  approach 

Lupas  and  Verdii  used  the  inverse  of  the  matrix  P  to  separate  the  signals  [2]. 
That  is; 

z  =  P~^x  =  Ah+P~^n  (4) 

=  Ab  +  ^  (5) 

where  z  =  [zi,  •  ■  zjf],  and  without  loss  of  generality  we  dropped  the  depen¬ 
dency  on  z.  It  follows  that  E  Apart  from  noise,  note  that 

z  contains  data  from  only  one  user.  Using  a  decision  on  the  elements  of  z, 
we  get  as  output  b  =  sgn(z).  Clearly  this  detector  is  near-far  resistant.  The 
probability  of  error  for  bk  depends  on  the  signal-to-noise  ratio 


SNRk  = 


ai  _  ai  |Fl 


(6) 


where  |F|  is  the  determinant  of  F  and  lFfc|  is  the  co-factor  of  F. 

The  Verdii  separator  is  only  one  of  many  possible  separators.  In  fact  any 
linear  transformation  that  diagonalises  matrix  F  will  do.  In  the  next  section  we 
present  such  a  transformation  and  examine  its  properties  and  the  performance 
of  the  corresponding  multi-user  detector. 


Figure  1:  Transformation  of  matched  filter  outputs. 


3.2  Different  Decorrelating  Detector 


z  =  Vx  =VPAh+Vn 

=  VPAb+C 

We  will  try  to  find  V  ^  P~^  so  that  VP  is  a  diagonal  matrix.  We  propose 
to  use  V  =  I  —  W,  where 

0  u;i2  tuis  . . .  wix 

^21  0  W23  ■  ■  ■  W2K 

=  \  /ox 


d  WK-1,K 

L  .  .  •  WK,K-1  0 

as  shown  in  Fig.  1 

Note  that  W  is  not  necessarily  a  symmetric  matrix.  The  output  of  the 
detector  is: 

Zk=Xk-wlxk  +  Ck  (9) 

where  Wk  is  the  column  of  W  without  Wkk  and  Xk  is  the  vector  x  with  Xk 
taken  out.  We  propose  to  choose  W  so  that 

E[zkbk]  =  Q 

where  the  k^^  element  is  taken  out  of  b  to  form  6*. 

If  the  signal-to-noise-plus-interference  ratio  for  any  of  the  Zj,j  F  is  large, 
then  £'[2:j;Sgn(zj)]  =  E  [zkbj  (l  —  Pe^)].  Under  this  condition  (10)  is  an  ap¬ 
proximation  of 

E[zkSgn{zk)]  =  Q  ior  k=l,2,...,K  (H) 
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where  Zk  is  z  without  the  element.  By  using  (9)  we  get  for  (10): 

E  [(xk  - 'U)'[xk)bi!]  =  0  {oT  k  =  1,2, . .  .,K  (12) 

One  can  show  that 

Xk  =  Pk-Akbk  +  y/^hkPk  +  nk  (13) 


where  the  index  k  for  a  matrix  or  vector  means  that  rows  and  columns  k  are 
taken  out  of  the  matrix,  or  that  element  k  is  taken  out  of  the  vector.  Note  that 


E  [Kfchj:]  =  E  [{-^/^bk  +  PkAkbk  +  Tik)  bk] 

=  AkPk  (14) 

where  we  use  the  fact  that  the  data  of  the  different  users  are  uncorrelated  and 
independent  of  the  noise.  Also  using  (13)  we  have 


E  [bkxl 


E 


bk  {blAkPk  +  y/^bkPk  +nl^ 


AkPk 


(15) 


Therefore  substituting  (14)  and  (15)  in  (12)  we  get 

Wk  =  Pk^Pk  (^®) 

Using  these  Wk  the  linear  transformation  V  will  diagonalise  the  correlation 
matrix  P.  From  the  definition  of  W,  Pi  and  V  =  I  —  W 


VP  -  diag  (l  -  tnf  •  •  • ,  1  -  w'^pjn) 


and  for  the  noise  covariance 

E  [CC^]  =  E  [Vnn'^V'^]  =  a^VPV'^ 

VP  is  diagonal  and  the  diagonal  terms  of  V  are  unity,  hence 
E[Ck]  =  (i-'^kPk)<^^ 


Therefore, 

SNRk 


(1  -  Pk) 

«»■  (l  -  Pk  Pk) 


(17) 

(18) 


(19) 


(20) 


where  we  also  used  (16). 

Using  E  [Cj]  >  0  and  the  fact  that  Pk  is  positive  definite,  we  can  show  that 

0  <  1  -  Pfc  (-Pit)"  Vife  <  1  (21) 
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sgn(  z,) 


sgn(z,) 

Figure  2:  Decorrelator  for  two  users. 


Furthermore,  it  can  be  shown  that 

1  -  pI  {PkT^ 

Now  comparing  (20)  with  (6)  we  conclude  that  if  we  use  E[zkhk]  =  0  as 
a  diagonalising  transformation,  then  the  SNR  at  the  output  (and  hence  also 
the  error  probability  at  the  decorrelator  output)  is  equal  to  that  of  the  Verdii 
transformation. 

3.3  An  Improved  Decorrelating  Detector 

Using  (11)  directly  will  result  in  a  linear  transformation  that  will  not  separate 
the  signals  totally,  but  will  leave  some  interference  residue  in  the  output.  The 
additive  noise,  however,  will  be  smaller  so  that  the  output  signal-to-noise-plus- 
interference  ratio  will  be  better  (or  at  least  as  good)  than  the  total  separation 
signal-to-noise  ratio.  This  gives  us  improved  performance  in  the  low  interference 
region.  Solving  (11)  analytically  without  using  the  high  signal-to-noise  assump¬ 
tion  is  rather  difficult  due  to  the  correlated  noise.  However,  using  numerical 
solutions  we  find  that  particularly  in  the  region  where  interference  power  is  low, 
the  output  SNR  is  better  and  so  is  the  probability  of  error.  To  show  this  without 
relying  solely  on  numerical  methods  we  use  the  two-user  case.  Extending  the 
proof  to  a  multi-user  case  is  possible  following  similar  reasoning.  From  Fig.  2: 

=  \/^^l  + 

Xi  =  +  P-v/oTfei  +  n2  (23) 

with  E  [rij]  =  E  [n^]  —  and  £'[nin2]  =  pcr^.  Also 

Zi  —  Xi  —  11)12X2 

.  =  y/a^il-wi2p)bi+y/a^{p-wi2)b2  + 

ni  -  Wi2n2  (24) 
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and 


zi  —  ■  Xi  -  W21X1 

=  {I  -  W21p)  b2  +  {p  -  W21)  bi  + 

n2  -  W2ini  (25) 

The  variance  of  the  noise  in  output  zi  is: 

(^loise  Z,  =  -  2u;i2P  +  wl^)  (26) 

and  the  signal-to-noise  plus  interference  at  this  output  is: 

SNIR,,  =  — - (1  -^r^P)  - ^  (27) 

(T^  (1  -  2wi2P  +  ^12)  +  02  (p  -  W'12) 

If  oi  >  02  and  oi  >  then  SNIRz^  >  1  provided  that  1012  is  of  the  same 

order  as  p.  Under  these  conditions  sgn(zi)  ~  fei,  and  hence  £■  [22Sgn (zi)]  =  0 

will  lead  to  1021  =  p.  When  these  conditions  are  met,  let 

wy2  =  P  +  b  ,(28) 

To  find  8  we  must  equate 

E[zisgn{z2)\ 

=  E  [(-v/di  (1  -  p^  -  p8)  bi  -  8^fa^b2  +  C[)  ' 

(sgn(yh;(l-p'*)62  +  C2))]  (29) 

to  zero,  where  we  used  (24)  and  (25)  together  with  1021  =  P  and  u;i2  -  p  +  6, 
and  we  also  defined 

(^[  =  m  —  Wi2n2  —ni-{p  +  8)n2 

C2  =  02  -  tn2iOi  =712-  poi  (30) 

(29)  can  be  separated  into 

E  [(Voi  (1  -  p^  -  p8)  bi  -  by/a^b-^  ■ 

sgn  (yii'  (1  -  p^)  62  +  C2)]  + 

E  [C(sgn  (y^  (1  -  p^)  62  +  C2)]  (31) 

By  performing  the  expectation  it  can  be  shown  that 


1  +  ^fLSNRjj^  (1  -  2Q  (yiSNR)) 

where  LSNR  is  the  limiting-signal-to-noise  ratio  at  the  ou,tput: 

02  (1  —  p^) 

LSNR =  ^ 
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From  (32)  we  notice  that  6  has  the  opposite  sign  of  p  and  that  its  absolute 
value  |(5|  decreases  from  |p|  to  zero  when  LSNR  increases  from  zero  to  infinity. 

We  now  calculate  the  error  probability  at  zi.  Using  Wi2  given  by  (28),  and 
(24)  we  have 

zi  =  1/57  (1  ~  “  P^)  bi  +  y/^6b2  +  Cl 

where  Cl  is  given  by  (30) .  The  probability  of  error  can  be  found  to  be 


ei 


/  ysr  (1  -  p^ + \ 
2  Y  (7a/i  -  +  s^  j 


-Q 

2^ 


(T^/^  -  /)2  -1-  ^2 

\  (r\/l  -  -f  ^2  j 


+ 


(35) 


For  small  a2,  S  —p.  Then  Pei  Q  the  error  corresponding  to  a 

simple  BPSK  signal.  When  02  is  large,  6  — >  0,  and  Pej  — +  Q  (^y/a^\/l  —  p^/u’^ , 
the  error  corresponding  to  Verdu’s  decorrelating  detector. ' 


3.4  An  Adaptive  Version  of  the  Improved  Decorrelating 
Detector 

The  transformation  which  corresponds  to  the  solution  of  (11)  can  be  obtained 
as  a  steady  state  of  the  following  adaptive  algorithm,  previously  referred  to  in 
the  literature  as  the  ’’bootstrap  algorithm”  [10]: 

Wk{i+l)  =  Wk{i)+pzkSgn(zk)  (36) 

for  A:  =  1,2,..., P: 

With  this  algorithm  the  steady  state  output  error  probability  is  upper- 
bounded  by  a  value  equal  to  that  obtained  by  Verdu’s  decorrelator.  This  value 
is  reached  only  when  the  interference-to-noise  ratios  for  all  the  other  users  are 
large  in  comparison  to  the  desired-signal-to-noise  ratio. 


4  Numerical  Results 

Although  only  the  derivation  for  the  two-user  case  is  presented,  a  similar  analysis 
can  be  carried  out  for  the  multi-user  case.  Results  for  both  the  two-user  and 
the  multi-user  cases  are  presented. 

For  two  users.  Fig.  3  shows  the  dependence  of  (5  on  LSNR  for  different  values 
of  p.  The  corresponding  probability  of  error  is  given  in  Fig.  4.  It  can  be  seen 
from  these  graphs  that  the  improved  decorrelator  acts  at  least  as  good  as  the 
total  signal  separator,  going  asymptotically  to  the  results  of  the  separator,  for 
high  interference  power. 
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Results  for  multi-user  systems  are  shown  in  figs.  5  to  7.  For  a  system  with 
high  cross-correlation,  Fig.  5  shows  6n  and  613.  The  error  probability  for  user 
1  in  this  system  is  given  in  Fig.  6.  It  varies  between  the  error  probability  for  a 
single  user  BPSK  system  for  low  interference  and  Verdii’s  error  probability  for 
high  interference. 

A  simulation  was  run  for  a  system  with  cross-correlation  -1/7  for  all  users. 
In  Fig.  7  the  results  are  compared  to  the  theoretical  curve,  which  confirms  the 
outcome.  Again,  the  total  signal  separator  error  and  the  BPSK  error  are  shown 
to  bound  the  improved  decorrelator. 


Figure  3;  (5  as  a  function  of  LSNR  for  different  values  of  p. 
Eb/NOi  =  MB. 


Figure  4:  Theoretical  error  probability  of  user  1  as  a  func¬ 
tion  of  the  energy  of  user  2.  Eb/NOi  —  8dB. 
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EbA^O  users  2  &  3,  dB 


References 

[1]  S.  Verdii,  ’’Minimum  probability  of  error  for  asynchronous  Gaussian  mul¬ 
tiple  access  channels,”  IEEE  Transactions  on  Information  Theory,  Vol. 
IT-32,  No.  1,  Jan.  1986,  pp.  85-96. 

[2]  R.  Lupas  and  S.  Verdii,  ’’Linear  multiuser  detector  for  synchronous  code 
division  multiple  access  channels,”  IEEE  Transactions  on  Information  The¬ 
ory,  Vol.  IT-35,  No.  1,  Jan.  1989,  pp.  123-136. 

[3]  R.  Lupas  and  S.  Verdii,  ’’Near-far  resistance  of  multiuser  detectors  m  asyn¬ 
chronous  channels,”  IEEE  Transactions  on  Communications,  Vol.  COM-38, 
No.  4,  Apr.  1990,  pp.  496-508. 

[4]  M.K.  Varanasi  and  B.  Aazhang,  ’’Multi  stage  detector  in  asynchronous  code 
division  multiple  access  communications,”  IEEE  Transactions  on  Commu¬ 
nications,  Vol.  COM-38,  No.  4,  Apr.  1990,  pp.  509-519. 

[5]  M.K.  Varanasi  and  B.  Aazhang,  ’’Near-optimum  detector  in  synchronous 
code  division  multiple-access  systems,”  IEEE  Transactions  on  Communi¬ 
cations,  Vol.  COM-39,  No.  5,  May  1991,  pp.  725-736. 

[6]  Z.  Siveski,  Y.  Bar-Ness  and  D.W.  Chen,  ’’Error  performance  of  a  syn¬ 
chronous  multiuser  code  division  multiple  access  detector  with  a  multi 
dimensional  adaptive  canceler,”  The  European  Transactions  on  Telecom¬ 
munications  and  Related  Teehnologies,  Vol.  5,  No.  6,  Nov.-Dee.  1994. 

[7]  D.W.  Chen,  Z.  Siveski  and  Y.  Bar-Ness,  ’’Synchronous  multiuser  CDMA 
detector  with  a  soft  decision  adaptive  canceler,”  The  28‘’^  Annual  Confer¬ 
ence  on  Information  Sciences  and  Systems,  Princeton,  NJ,  Mar.  1994. 

[8]  Y  Bar-Ness  and  N.  Sezgin,  ’’Adaptive  Threshold  Setting  for  Multi-User 
CDMA  Signal  Separators  with  Soft  Tentative  Decision,”  The  29^''  Annual 
Conference  on  Information  Sciences  and  Systems,  Baltimore,  MD,  Mar. 
1995. 

[9]  S.  Verdii,  ’’Recent  progress  in  multiuser  detection,”  Advances  in  Commu¬ 
nication  and  Signal  Processing,  Springer  Verlag,  1989. 

[10]  A.  Dine  and  Y.  Bar-Ness,  ’’Bootstrap:  a  fast  unsupervised  learning  al¬ 
gorithm,”  IEEE  1992  International  Conference  on  Acoustics,  Speech  and 
Signal  Processing,  ICASSP  ’92,  San  Francisco,  CA,  Mar.  23  1992,  paper 
43.8. 


41 


Appendix  B 


Comparative  Study  of  the  Linear  Minimum  Mean 
Squared  Error  (LMMSE)  and  the  Adaptive 
Bootstrap  Multiuser  Detectors  for  CDMA 
Communications  * 


Hongya  Ge  and  Yeheskel  Bar-Ness 
Center  of  Communications  and  Signal  Processing  Research, 
Dept,  of  Elec.  Sz  Comp.  Engg. 

New  Jersey  Institute  of  Technology, 

University  Heights,  Newark,  NJ  07102 


Abstract 

In  this  work,  we  analyze  the  performance  of  linear  minimum  mean 
squared  error  (LMMSE)  estimate-based  multiuser  detector  for  CDMA 
Communication  Systems.  Through  numerical  evaluation  and  computer 
simulations,  we  also  compare  the  performance  of  the  LMMSE  detector 
[1]  and  the  adaptive  bootstrap  multiuser  detector  [2].  Our  analysis  and 
simulations  show  that  even  though  these  two  detectors  were  proposed 
based  on  different  optimization  criteria,  they  exhibit  approximately  equal 
performance  in  multiuser  CDMA  communications  applications. 


1  Introduction 

Multiuser  separation  and  interference  suppression  is  an  active  research  topic  in 
CDMA  communications  area.  Various  detectors  have  been  proposed  to  balance 
the  computational  simplicity  and  reliable  detection  performance  [1  —  7].  Most 
of  the  proposed  detectors  treat  the  multiuser  signal  vector  as  a  deterministic 
vector,  and  work  only  on  the  conditional  probability  density  function  (PDF)  of 
the  data,  given  the  multiuser  signal  vector.  We  notice  that  the  information  bit 
of  multiuser  is  actually  a  random  vector.  Therefore  the  optimal  approach  of 
separating  the  information  bits  of  multiuser  and  suppressing  the  interference  is 
equivalent  to  estimating  each  realization  of  the  multiuser  information  bits  with 

‘This  work  was  supported  in  part  by  the  Office  of  Sponsored  Research,  NJIT,  and  Rome 
Air  Force  Lab  under  contract  F30602-94-C-0135 
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minimum  probability  of  error.  By  incorporating  a  prior  knowledge  (in  a  statis¬ 
tical  sense)  of  the  multiuser  signal  vector,  one  can  always  improve  the  overall 
detection  performance.  The  availability  of  the  a  prior  knowledge  depends  on  the 
specific  communication  systems.  Under  the  linearity  constraint,  we  investigated 
the  linear  minimum  mean  squared  error  (LMMSE)  estimate-based  multiuser  de¬ 
tector,  analyzed  its  detection  performance,  and  compared  with  the  performance 
of  the  adaptive  bootstrap  multiuser  detector  [2] .  It  should'  be  pointed  out  that 
the  LMMSE  detector  was  proposed  in  [1]  and  was  termed  as  MMSE  detector. 
Based  on  statistical  signal  processing  concept,  the  MMSE  estimate  is  obtained 
by  finding  the  conditional  mean  of  the  posterior  PDF  of  the  multiuser  infor¬ 
mation  bit,  given  the  received  data.  And  this  MMSE  estimate-based  detector 
is,  in  general,  a  nonlinear  detector.  If  we  constrain  our  detector  to  be  in  the 
linear  class,  we  can  find,  as  proposed  in  [1]  and  not  accurately  termed  as  MMSE 
detector,  the  sub-optimal  linear  MMSE  (LMMSE)  estimate-based  detector.  In 
this  paper,  we  further  analyze  the  performance  of  the  LMMSE  detector  in  detail 
and  verify  our  results  through  numerical  evaluation  and  computer  simulation. 
Comparative  study  of  performances  of  the  linear  class  of  decorrelating  detector, 
LMMSE  detector,  and  adaptive  bootstrap  multiuser  detector  is  provided  in  this 
work. 


2  Description  of  the  Problem 

Due  to  the  multiple  access  (MA)  scheme  used  in  CDMA  systems,  the  avail¬ 
able  baseband  data  r{t)  at  the  receiver  is  actually  a  mixture  of  multiuser  data 
embedded  in  additive  noise.  That  is, 

K  _ 

m  ^'t(*)  Sk{t  -iT-Tk)  +  n{t) ,  (1) 

i  k  =  l 

where  K  is  the  number  of  users,  ak{i),  bk{i),  Sk{t),  and  Tk  are  the  bit  energy, 
information  bit,  signature  waveform,  and  transmission  delay  of  the  kih  user  in 
the  ith  symbol  interval  (of  duration  T),  respectively.  n{t)  is  a  white  Gaussian 
process,  with  two-sided  power  spectral  density  of  cr^. 

In  this  work,  we  only  consider  the  case  when  r*,  =  0  (fc  =  1, . . . ,  K),  which 
corresponds  to  the  synchronous  channel.  The  synchronous  channel  model  is 
valid  for  down-link  channel  (from  base  station  to  mobiles).  Note  that  once  the 
channel  is  synchronized,  all  the  information  bits  of  the  multiusers  in  the  ith 
symbol  interval  are  completely  contained  in  the  data  r{t)  within  the  ith  symbol 
interval.  Therefore,  we  can  concentrate  on  solving  multiuser  separation  problem 
within  a  specific  symbol  interval  (say  i  =  0  th  interval)  without  lose  of  generality. 

Let  us  consider  a  specific  bit  interval,  say  i  =  0.  After  ignoring  index  i,  we 
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have, 

r{t)  =  [  si{t)  S2{t)  ■■■  s^it)  ]  ■  A-h  +  n{t) ,  0<t<T,  (2) 

■  6i 
^2 

with  A  =  diag{  ^/a^,  ■  •  • ,  and  b  = 

L 

We  then  filter  r(t)  with  a  bank  of  matched  filters,  whose  impulse  responses 
are  given  by, 

hk{t)  =  SkiT-t)  ,  k  =  l,2,...,K. 

Columnize  the  sampled  outputs  of  the  bank  of  matched  filters  at  t  =  T,  we  get 
the  following  matrix  notation, 

X  =  P  •  A  •  b  +  n ,  (3) 

where  x  =  [  a;i(T)  X2{T)  ■  •  ■  Xk{T)  ]^,  with  Xk{T)  —  r{t)  *  hk{t)  \t=T  be¬ 
ing  the  feth  matched  filter  output  sampled  at  time  instant  t  —  T\  and  n  ~ 
<T^  P)  being  the  colored  Gaussian  noise  due  to  the  matched  filtering.  Note 
that  the  model  in  (3)  is  also  valid  for  asynchronous  channel,  except  for  a  larger 
dimension.  Since  for  asynchronous  channel,  by  introducing  the  partioned  signa¬ 
ture  waveforms  si(t),  sf  (0)  ' ' '  >  ^nd  s^{t)  as  in  [3],  the  resultant 

matrix  P  in  (3)  is  of  a  dimension  (2K  —  1)  x  (2K  —  1).  If  the  matrix  P  is  nonsin¬ 
gular,  one  may  use  the  one-shot  method  of  [3].  For  more  general  asynchronous 
channel,  we  also  proposed  a  multi-shot  approach  in  [7]  to  handle  the  possible 
singularity  problem.  In  the  sequel,  we  only  consider  the  synchronous  channel 
for  notational  simplicity. 

Note  that  matrix  P  in  (3)  is  the  cross-correlation  matrix  of  the  signature  wave¬ 
forms.  P  is  symmetric,  positive  definite  and  nonsingular,  and  its  elements  are 
given  by, 

P[i,  j]  =  /  Si{t)  Sj{t)  dt  =  pij  .  \ 

In  practice,  due  to  the  finite  bandwidth  constraint  and  the  existence  of  a  large 
number  of  users,  the  signature  waveforms  are  not  idealy  orthonormal.  Therefore 
the  matrix  P  will  not  be  an  identity  matrix  in  general.  The  non-diagonal  nature 
of  the  P  matrix  will  cause  the  multiple  access  interference  (MAI).  In  order  to 
remove  the  MAI,  various  detectors  have  been  proposed  [1  —  7].  One  major  effort 
of  proposing  various  detectors  is  trying  to  balance  the  computational  simplicity 
and  reliable  detection  performance.  In  the  following  sections,  we  first  summarize 
various  linear  multiuser  detectors,  and  propose  new  approa|Ches  to  analyze  their 
performances  from  statistical  signal  processing  and  parameter  estimation  point 
of  view,  we  then  compare  their  performances  through  numerical  evaluations 
and  computer  simulations. 
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3  Linear  Class  of  Multiuser  Detectors 

In  this  section,  we  review  some  recently  proposed  linear  multiuser  detectors 
from  statistical  signal  processing  point  of  view,  analyze  their  performances,  and 
make  a  comparative  study  of  the  decorrelating  detector,  the  LMMSE  detector, 
and  the  adaptive  bootstrap  multiuser  detector. 


3.1  Simple  Decorrelating  Detector 

The  simple  decorrelating  detector  [5]  is  originated  by  finding  a  linear  condi¬ 
tional  maximum  likelihood  estimate  (MLE)  of  the  signal  vector  0  =  A  •  b  from 
the  conditional  PDF  p(x|e)  obtained  from  (3).  It  then  detects  the  multiuser 
information  bit  b  by  directly  making  decision  on  the  linear  conditional  MLE 
B  =  •  X,  which  is  simply  a  linear  matrix  inverse  filtering  operation  on  the 

data  X.  The  outputs  of  the  multiuser  decorrelator  can  be  written  as, 

z  =  P”' -x  =  A -b -h  V,  (4) 

with  V  ~  ^”(0,  <T^P"‘)  being  the  colored  Gaussian  noise  with  a  different  co- 
variance  matrix. 


From  (4)  we  can  see  that  the  signal  bits  are  separated.  The  detector  then 
makes  decision  from, 

b  =  sign(z) .  (5) 

This  detector  has  the  advantage  of  structural  simplicity.  It  is  also  near-far 
resistant.  But  a  potential  problem  associated  with  this  detector  is  that  noise 
is  enhanced  by  the  P“'  inverse  filtering.  We  show  that  this  detector  has  its 
limited  performance  from  the  following  simple  analysis.  Content-Length:  19623 
X-Lines;  461  Status;  O 

Assume  that  we  are  interested  in  getting  the  probability  of  error  expression  of 
the  itth  user  of  the  decorrelating  detector  in  (5),  where  vector  z  is  defined  in  (4) 
and  Ab  is  assumed  fixed.  We  can  then  re-arrange  the  order  of  the  components 
of  data  vector  z  in  (4)  as  following, 


A 

Zk 

^k 

4- 

Vk 

Z  — 

.  -k  . 

Afcbjt 

.  -k  . 

where  bj,  is  a  (K  —  1)  dimensional  vector  constructured  from  the  original  K 
dimensional  vector  b  with  its  fcth  component  removed;  Ak  is  a  {K  - 1)  x  (A  - 1) 
matrix  constructured  from  the  original  K  x  K  dimensional  matrix  A  with  row 
and  column  associated  with  the  feth  user  removed.  Similar  notations  also  apply 
to  Zj.,  Y-ki  Pk’  sections.  When  Ab  is  given  and  fixed, 

we  have. 


Zk 

U 


~  h! 


^/^k 

Afcbjt 


(7) 
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Note  the  fact  that 


1  Bkl] 
Pk  P* 

—1 

1 

i-i'rp. 

Hence,  for  a  given  A  b,  ~  A/" 
of  the  A:th  user  can  be  easily  found  as, 


probability  of  error 


Pe{k)  =  P{zk>0\bk  =  -1}  P{bk  =  -1}  + 
P  {^k  <  0  I  =  — h}  P  {bk  =  +1} 


(8) 


where  Q(x)  =  /  -y=  e  *  dt,  and  we  assume  that 

Jr 

P  {bk  =  +1}  =  P  {bk  =  -1}  =  1/2 . 

Note  that  in  (8),  the  matrixP*  is  a  (A— l)x(A— 1)  matrix,  constructed  from 
P  matrix  after  removing  the  contribution  of  the  fcth  user.  Matrix  Pj,  is  positive 
definite,  since  it  is  a  sub-matrix  of  the  positive  definite  matrix  P.  Therefore 
the  inequality  0  <  (1  —  Pj *  Pj.)  <  1  will  always  hold,  with  the  equality  holds 
if  and  only  if  P  is  a  diagonal  matrix.  This  corresponds  to  the  case  of  using  a 
set  of  perfect  orthonormal  signature  waveforms.  Hence,  the  performance  of  the 
decorrelating  detector  in  (5)  will  always  worse  than  the  BPSK  limit,  as  shown 
in  the  inequality  in  (8).  The  near-far  resistance  of  the  decorrelating  detector 
can  be  easily  seen  from  its  Pe{k)  expression  in  (8),  since  Pe{k)  is  invariant  to 
Ajfc. 

We  further  notice  the  following  fact  that  the  above  decorrelating  detector 
is  also  a  linear  estimate-based  detector.  Part  of  its  limited  performance  is  due 
to  the  fact  that  linear  estimate  6  =  P~’  x  is  based  on  the  conditional  PDF 
p(x|0).  Therefore,  we  can  expect  to  further  improve  its  detection  performance 
and  maintain  its  linear  feature,  by  incorporating  the  joint  statistics  of  both  6 
and  X  into  the  estimate  cis  demonstrated  in  the  following  linear  minimum  mean 
squared  error  (LMMSE)  estimate-based  detector. 
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3.2  Linear  Minimum  Mean  Square  Error  (LMMSE)  Estimate- 
based  Multiuser  Detector 

It  is  well  known  that  among  the  linear  class  of  estimates,  LMMSE  estimate 
exhibits  the  minimum  mean  squared  estimation  error .  Therefore  we  can  improve 
the  detection  performance  of  the  above  linear  decorrelating  detector  by  deriving 
a  LMMSE  estimate-based  multiuser  detector.  Specifically,  let  us  rewrite  the 
matched  filter  outputs  in  (3)  as, 

x=:PAb-bn  =  P0-t-n.  (9) 


We  assume  that  all  the  components  of  the  random  multiuser  information  bit 
b  are  independent  and  identically  distributed  (i.i.d.)  with  zero  mean  and  unit 
variance.  Even  further,  the  random  vector  6  and  noise  vector  n  are  statisti¬ 
cally  independent.  For  most  communication  applications,  these  assumptions 
are  reasonable  ones.  We  then  define  a  new  random  vector  as  follows. 


A 

■  9  ■ 

■  Ab  ■ 

A  o' 

b  ' 

y  = 

X 

— 

X 

PA  I 

n 

s 

✓ 

H 


The  expectation  and  covariance  matrix  of  the  new  vector  y  can  be  found  as, 

My)  =  o, 


cov(y)  =  H  cov( 


b 

n 


)H’ 


A^P 

P  A^  P  A^  P  +  P 


L  Y^-x.6  Sxx  J 

Then  LMMSE  estimate  of  9  can  be  formed  from  the  following  [8,9], 


9 LMMSE  —  ^{^)  +  120X  ^xx 
=  1^0  X  I^XX 

=  (P-1-<t2A-2)“^  x. 

' - - - ^ 

W 

The  LMMSE  estimate-bsaed  multiuser  detector  makes  its  decision  based  on  the 
decision  rule, 

b  =  sign  (  9lmmse  )  •  (1^) 
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Note  that  the  diagonal  matrix  <7^  A  ^  involved  in  W  of  (10)  is  actually  the 
inverse  SNR  matrix^  i.  e. 


(T^  A 


-2 


ai/a^  0 
0  02/0-2 


0 

0 


-1 


0  0  •  •  •  aKjiP' 


(12) 


When  the  interferences  of  other  users  are  small  compared  to  the  noise  level,  W 
reduces  into  a  diagonal  matrix,  or  W  =  (  P  +  o-^  A“2)-i  r;  (I -f  A“2)-i. 
In  this  case,  LMMSE  detector  has  the  same  performance  as  that  of  the  single 
user  detector,  which  is  the  BPSK  limit.  When  the  interference  levels  are  very 
large  compared  to  the  noise  level,  then  W  =  ( P  +  o-^  A~2  )-i  ^  p-i^  LMMSE 
detector’s  performance  is  compariable  to  that  of  the  decorrelating  detector. 
Therefore,  the  overall  performance  of  LMMSE  detector  is  better  than  that  of 
the  decorrelating  detector.  Even  further,  we  can  calculate  the  probability  of 
error  of  the  I;th  user  based  on  the  following  observation. 

We  arrange  the  order  of  all  the  users  such  that  9'^  =[  6^  ]  •  We  then 

decompose  the  above  derived  Bimmse  as  follows. 


^LMMSE  —  Wx=:  (P  +  X 

=  (P<t2a-2)“^  •  (PAb  +  n) 
=  6>-<72WA-ib  + Wn 


(13) 


=  9  +  e 


where  the  estimation  error,  e  —  — cr^  W  A  ^  b  +  W  n,  contains  both  bias  and 
noise  components.  We  will  notice  later  that  the  improved  performance  of  the 
LMMSE  detector  is  achieved  by  trade-off  a  little  bias  for  less  noise  variance, 
which  finally  results  in  less  overall  mean  squared  error  (M^E). 

For  a  given  information  bit  b,  LMMSE  detector  of  (10)  makes  an  errorous 
decision  on  the  I;th  user’s  information  bit  bk  whenever, 

OR  ^  =  +1, 

e*  >  +y/^;  when  6;;,  =  -1- 
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Therefore,  the  probability  of  error  for  the  kth  user  can  be  expressed  as, 

P,{k)=  Peik\b,,h^)  P{bk)Pih,) 

=  P{ek>+^/b^\bk  =  -hhc)P(h  =  -l)+ 

Pick  <  I  bk  =  +1,  b,)  Pibk  =  +1)  P  {h^} 

=  ^  {P{ek  > +\^\bk  = -l,hk)  + 

P{ek  <  \bk  =  +1,  bj.)  }  •  P  (bjt) 

(14) 

where  Rk{-)  =  A"^  b,  o-|  =  cr^  w[  P  w*,  and  v/l  is  the  ^th  row  of  W 

matrix.  Or, 

[1  -pl{Pk  +  <r^A-,Y^] 

1  +  o'^/at  -  P^(Pit  +  pj. 

We  also  verified  the  fact  that  under  weak  and  strong  interference  conditions,  the 
above  error  probability  expression  will  reduce  either  into  the  single  user  BPSK 
limit  or  into  that  of  the  decorrelating  detector  as  follows. 


Peik)=< 


/akil-pTPl^P,) 


small  interferences 


,  strong  interferences 


For  the  case  of  K  —  2,  the  above  expression  can  be  simplified  as. 


iZl(-l,+l)  =  (T' 


i2l(-l,-l)  =  <7 


2  (l  +  <T^/a2)  /v^  +  p/y/^ 

(1  +  <72/ai)(l  +  <72/02)  -1/92 

2  (1  +  <7^/02)  /yaT  -  p/y/^ 
(1  +  (72/ai)(l  +  <72/02)  -  /9^ 


<7i  =  (7  ■ 


'((1  +  <72/02)  -  p2)^  +  p2  (1  _  p2) 

(1  +  <72/ai)(l  +  <72/02)  - 
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We  can  also  easily  verify  that  the  following  facts, 

lim  Pe(l)  =  q{—\  ,  BPSK  limit 

^  ^  ^  ^  ^  ^  detector 

3.3  Adaptive  Bootstrap  Multiuser  Detector 

In  [2],  an  adaptive  bootstrap  multiuser  detector  was  proposed  by  minimizing 
the  correlation  between  the  canceler  output  of  the  desired  pser,  say  Zk,  and  the 
hard-limiter  decision  on  the  canceler  output  of  other  users,  say  sign(zj(,).  That 
is,  the  matrix  W^,  with  its  diagonal  elements  being  zeros,  is  adaptively  chosen 
such  that 

£'{2i:sign(zjfe)}  =  0,  1,  2,...,^!:,  (16) 

r  "xT 

where  z  =  [  Zk  zj  ]  =  (I  —  W)  x .  The  adaptive  version  of  the  alogrithm 

updates  the  A;th  column  of  W  with  its  ^th  zero-element  removed,  say  Wt  as 
follows, 

i5Lfc(i+l)  =  w*(0  +  /iZtsign(z;fe),  k  =  l,2,...,K.  (17) 

Once  matrix  W  is  obtained  adaptively,  the  bootstrap  detector  makes  its  decision 
based  on  the  decision  rule 

b  =  sign  ((I  -  W)x)  . 

The  probability  of  error  of  the  kth  user,  using  an  adaptive  bootstrap  multiuser 
detector,  was  calculated  in  [2]  as, 

^  ,  (18) 

b. 

with  6  =  + 

‘^\/i-£iprp,+£prp, 

4  Numerical  evaluations  and  computer  simula¬ 
tions 

Numerical  evaluation  and  computer  simulation  results  in  this  section  further 
confirm  our  above  derivations.  Note  that  implementation  o^the  LMMSE  detec¬ 
tor  needs  information  on  matrix  P  and  SNR  matrix  of  (11),  but  the  adaptive 
bootstrap  multiuser  detector  developed  by  BarNess  et.  al.  [2]  can  achieve  the 
same  performance  as  that  of  the  LMMSE  detector  without  these  information. 
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In  the  following  figures,  we  show  the  performance  comparison  of  the  decorre- 
lating  detector,  the  LMMSE  detector,  and  the  adaptive  bootstrap  multiuser 
detector  for  high,  medium,  and  low  capacity  cases.  We  also  plot  the  BPSK 
limit  as  a  reference  lower  bound  on  the  performance.  It  can  be  seen  that  there 
is  an  equivalence  between  the  LMMSE  detector  and  the  adaptive  bootstrap 
multiuser  detector.  The  improved  performances  of  both  the  LMMSE  and  the 
adaptive  bootstrap  multiuser  detectors  are  obtained  by  the  idea  of  trade-off  a 
little  bias  for  less  noise  variance,  which  finally  results  in  an  overall  less  mean 
squared  error  (MSE).  And  the  adaptive  bootstrap  multiuser  detector  provides 
a  practical  implementation  of  the  LMMSE  detector. 


5  Conclusions 

We  investigate  the  performance  of  the  /inear  minimum  mean  squared  error  (LMMSE) 
estimate-based  detector  and  the  adaptive  bootstrap  multiuser  detector.  Nu¬ 
merical  evaluation  and  computer  simulations  verify  our  derailed  analysis  on  the 
performance  of  the  LMMSE  detector.  The  theoretic  analysis  and  computer  sim¬ 
ulations  also  confirm  the  inherent  relationship  between  the  LMMSE  and  adap¬ 
tive  bootstrap  detectors.  For  the  consideration  of  practical  implementation,  the 
adaptive  bootstrap  multiuser  detector,  which  does  not  need  the  knowledge  of 
the  multiusers’  signal-to-noise  (SNR)  matrix,  is  preferred  for  CDMA  communi¬ 
cations  applications. 
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Pe(11  versus  SNHfi).  (i=2. 3) 


(a).  High  capacity  case 


Figure  1;  Performance  of  the  LMMSE  detector  and  the  adaptive  bootstrap 
detector  in  multi-user  CDMA  system.  Also  shown  in  figure  are  the  results  of  the 
decorrelating  detector  and  the  BPSK  limit.  Parameters  used;  SNRi  =  8  dB, 
K  =  3,  M  =  1000000  independent  trials. 
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Appendix  C 


Adaptive  Multiuser  Bootstrapped  Decorrelating 
CDMA  Detector  for  One-shot  Asynchronous 
Unknown  Channels  * 

Yeheskel  Bar-Ness  Nadir  Sezgin 


Abstract 

Vast  research  was  recently  performed  on  signal  detection  of  multiuser 
Code  Division  Multiple  Access  (CDMA).  Particularly  for  uplink  (users  to 
base  station)  the  signals  are  asynchronous  and  the  near-far  problem  is  an 
important  issue  to  deal  with.  All  near-far  resistant  detectors,  adaptive 
or  non-adaptive,  assume  knowledge  of  the  relative  delays  of  the  differ¬ 
ent  users’  signals.  Among  these  are  the  one-shot  detectors  suggested  by 
Verdu.  In  this  paper  we  suggest  an  adaptive  algorithm  to  decorrelate  the 
outputs  of  the  one  shot  matched  filters  that  assumes  no  knowledge  of  the 
relative  delay.  The  performance  of  this  approach  is  shown  to  be  better 
than  the  non-adaptive  “zero-forcing”  method  previously  used  and  which 
implements  linear  transformation  via  the  inverse  of  matched  filters’  output 
cross-correlation.  For  simplicity  a  two-user  case  is  presented.  Extension 
to  a  higher  number  of  users  is  relatively  simple. 

I.  INTRODUCTION 

CDMA  is  considered  to  be  a  promising  multiplexing  method  for  multiuser  per¬ 
sonal,  mobile  and  indoor  communications.  One  of  the  problems  a  designer  of 
such  a  system  is  faced  with  is  the  so  called  “near-far  problem,”  resulting  from 
excessive  Multiple  Access  Interference  (MAI)  energy  from  nearby  users,  com¬ 
pared  with  the  desired  user’s  signal  energy.  Power  control,  that  is,  adjustment 
of  transmitter  power,  depending  on  its  location  and  the  signal  energies  of  the 
other  users,  has  been  suggested  as  a  solution  to  this  problem.  But  it  requires 
a  significant  reduction  in  the  signal  energies  of  the  strong  users  in  order  for 
the  weaker  users  to  achieve  reliable  communication.  This  results  in  an  overall 
reduction  in  communication  ranges. 

An  optimum  near-far  resistant  multiuser  detector,  without  a  need  for  power 
control  has  been  proposed  by  Verdu[l].  Its  complexity,  however,  is  exponen- 

This  work  was  partially  supported  by  a  greint  from  Rome  Air  Force  Lab  (AFSC,  Griffss 
Air  Force  Base)  imder  the  contract  F30602-94-C-013S 
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tial  in  terms  of  the  number  of  users,  which  makes  it  unsuitable  for  practi 
cal  situations.  Later,  a  sub-optimum  detector  whose  complexity  is  linear  was 
proposed[2].  However  this  detector,  like  Verdu’s  optimum  detector,  needs  the 
knowledge  of  user  signal  energies. 

Adaptive  near-far  resistant  CDMA  detectors  that  do  not  need  the  knowledge 
of  the  received  signal  energies  have  been  developed  at  NJIT  in  recent  years[3-6]. 
All  these  detectors  assume  a  non-fading  multiuser  environment. 

All  asynchronous  detectors  (except  [6]),  adaptive  or  non-adaptive,  assume 
knowledge  of  the  different  users  signals  relative  delays. 

In  this  paper  we  propose  to  use  the  decorrelation  algorithm  for  separation 
of  the  user  signals  obtained  at  the  outputs  of  Verdu’s  one-shot  bank  of  matched 
filters.  For  this  we  do  not  require  knowledge  of  the  relative  delays. 

In  the  next  section  we  present  the  system  model  and  drive  the  relation 
between  the  matched  filters’  output  and  the  users  data.  Since  the  delays  are 
not  known,  the  cross-coupling  matrix  is  considered  unknown.  The  decorrelating 
algorithm  and  the  decorrelator’s  steady-state  weights  and  outputs  are  given  in 
Section  III.  In  Section  IV.  we  discuss  performance  evaluation  and  in  Section  V. 
we  present  an  example  and  results. 

II.  SYSTEM  MODEL 

For  the  asynchronous  multiuser,  the  equivalent  low-pass  signal  at  the  input  of 
the  matched  filter  bank  is  given  by. 


K 

~  +  T  —  Tk)  +  n{t),  (1) 

fc  =  l  i 

where  K  is  the  number  of  users,  ak,0k,bk,  Sk  and  Tk  are  the  signal  amplitude, 
carrier  phase  shift,  user  bit,  signature  waveform  and  the  relative  delay  of  the 
ifcth  user.  n{t)  is  the  zero  mean  AWGN,  with  a  two-sided  power  spectral  density 
ofiVo/2. 

For  the  sake  of  simplicity  we  will  restrict  ourselves  to  two  users  only.  Ex¬ 
tension  to  a  higher  number  is  relatively  simple. 

Representing  the  signal  of  Eq.  (1)  in  one-shot  of  the  ith  bit  of  user  one,  and 
without  loss  of  generality,  letting  i  =  0,  we  write, 

r(t)  =  +  ^/a2^-y=S2{t)b2{—l) 

+\/ ^2(1  “  f2)-"/==^=®2* (0^2(0)  +  ^(^)i  (2) 

y/l  —  (2 

where  0  <  r2  <  T  is  considered  unknown  or  estimated  with  an  error,  and; 


S^2it) 


S2it  +  T-r2)  iiQ<t<T2 
0  iiT2  <t  <T 
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0  if  0  <  <  <  r2 

S2{t-T2)  if  r2  <  f  <  T 


rT2 

(2—1  S2{t  +  T  —  T2)dt . 
Jo 


As  shown  in  Fig.  1,  we  apply  r(t)  to  a  bank  of  filters  by  using  as  a  second 
input  sf  (t),  respectively,  the  signature  function  of  Eq.  (3)  extended 

to  T2  instead  of  ,  normalized  by  €2  instead  of  €2  ■ 

The  output  of  the  first  filter, 

a:'i(0)  =  1/3761(0)  +  \/57e2/>i262(0)  +  n/O),  (4) 

where  ni(0)  =  n(t)si{t)dt  is  a  zero-mean  Gaussian  variable  with  variance  of 
No/2,  and 

1 

P21  =  -7=  /  Si(t)s2{t  +  T  —  T2)di 

v^2  Jo 

1 

~  yp,  J  -  r2)dt.  (5) 

The  subscript  in  x'(0)  indicates  that  we  are  using  zero  bit  of  user  one  in  our 
one-shot. 

The  output  of  the  second  filter  is  given  by. 


4(0)  =  /  r{t)s^'{t)dt 

v/f2  Jo 


rr  /  '  y’)^2 

Vf2  Jo 

\/^P21^l(0)  +  \/a2C2P22^62(  — 1) 
+  \/ a2(l  -  f2)P22^62(0)  -f  ^2(0), 


where 


1 

4i  =  —j=f  /  Sx{t)s2(t T  -  T2)dt 

V  (2  Jo 

L'L  1  . 

0^2  =  S2{t  +  T-T2)s2it  +  T-T^)dt 

L'R  1 

4"*  =  -J===  S2{t-T2)s2{t  +  T-T!2)dt 

V  ^2(1  ~  (2)  Jo 

and 

2(0)  =  /  n{t)s2{t-\-T-r!2)dt. 

V^2  Jo 
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The  output  of  the  third  filter  is  given  by, 


where 


Pi2  — 


.R'l  _ 
P22  — 


.R'R  _ 

P22  — 


n3(0)  = 


rp 

=  J—f  I  r{t)s^' {t)dt 

V  1  ^2  ‘'0 

+  1/02(1  —  €2)P22^^2i^)  +  ^3(0)1 

T 

-2==  f  Sl{t)s2it  -  T^)dt 

a  —  €2  Jo 

,  ^  _=  r  S2it  +  T-  T2)S2{t  -  T^)dt 

/(I  -  €9^2  Jo 

1 

.  /  S2(^  -  '’■2)S2(^  -  T2)Ji 

/(I  -  €2)(1  -  ^2)  Jo 

T 

1=  /  n{t)s2{t  -  'r^)d<. 

/I  -  e'o  -/o 


Combining,  Eqs.  (4),  (6)  and  (8)  in  matrix  form, 

x'(0)  =  P^Ab(0)  +  n(0), 

where 


1  P21  P12 

=  P2I  P22^  P22^  > 

.  P'l2  P22^  P22^  ■ 

A  =  diag[-/ar,  ^/d^,  y/ 02(1  —  £2)]  and 

b(0)  =  [6i(0),62(-l),fc2(0)]^  .  n(t)  is  zero-mean  Gaussian  vector  with  covari- 
ance, 

1  P2I  Pl2  W.  ^  ^ 

Rn=  P21  1  0  Y- 

.P12O  1  J  ^ 

2.1.  Known  Delay 

If  the  relative  delay  T2  is  known  then  we  choose  T2  =  T2,  and  as  a  result  P21  = 

P21,/l2  =  P12,P22^  =  P22^  =  P22^  =  P^^""  =  ^  4  =  £2  =  1-  The 

matrix  Pd  is  known.  If  |p2i|  +  IP12I  <  1  Pd  is  invertible.  A  simple  linear 
transformation  by  P^^  results  in, 

z(0)  =  Ab(0)  +  Po’’n(0),  (13) 

where  Pd  is  the  modified  cross-coupling  matrix.  Since  the  components  of  b(0) 
are  assumed  uncorrelated  then  the  components  of  z(0)  are  also  so. 


III.  THE  DECORRELATING  ALGORITHM 


3.1.  Unknown  Delay 


For  this  case  the  matrix  P  is  unknown.  We  propose  to  use  an  adaptive  algorithm 
to  decorrelate  the  outputs  of  bank  of  matched  filters  (see  Fig.  2).  From  this 
figure  we  have 


z'(0)  =  x'(0)  -  W^x'(O), 


(14) 


where 


W  = 


0  Wi2  Wiz 
W21  0  W2Z 

Wzi  WZ2  0 


The  kth  output  of  the  decorrelator  can  be  expressed  as, 


(15) 


Zk=Xk-'wlxk  (16) 

where  wj,  is  the  fcth  column  of  W  with  the  element  Wkk  deleted  and  x^  is  the 
vector  obtained  from  x  by  deleting  Xk- 

The  multidimensional  decorrelator  is  the  same  as  that  frequently  used  in 
neural  networks  and  other  applications  of  signal  separation[7-8]. 

For  controlling  the  weights,  we  use  the  steepest  descent  algorithm  which 
simultaneously  reduces  the  absolute  value  of  the  correlation  between  the  outputs 
of  the  decorrelator  and  the  decision  on  all  other  outputs.  That  is,  the  weight, 
Wik  is  controlled  by  the  recursion, 


Wmi  *-  Wmi  -  n{zisgn{zm)),  1  <  /,  m  <  3,  /  9^  m.  (17) 

Recursion  in  Eq.  (17),  reaches  steady-state  in  the  mean  when  £'{ziSgn(2:m)}  =  0. 
Note  that  W21,  for  example,  is  used  to  cancel  residue  of  62(— 1)  at  the  output  of 
zi.  It  will  settle  down  only  when  that  residue  (being  correlated  with  62(— 1))  is 
zero.  But  any  reduction  of  62(— 1)  at  zi  will  improve  6i(0)  (smaller  error)  and 
hence,  will  be  more  effective  in  reducing  the  residue  of  61  (0)  at  Z2  and  zz  through 
the  control  weights  W\2  and  Wiz-  Therefore,  the  process  of  residue  cancellation 
is  enhanced  successively  which  justify,  using  the  name  “bootstrap”  in  the  past. 
In  the  notation  of  Eq.  (16),  we  write  Eq.  (17)  in  vector  form, 

w*  <-  Wfc  -  //(2isgn(zjfe)),  (18) 

where  again  z*  is  obtained  from  the  column  of  z  by  deleting  Zk-  The  steady 
state  is  reached  if  E{zkSgn{zk)  =  0  for  F  =  1, 2, 3. 

Now,  sgn(zjfe)  =  bjfe,  and  if  we  assume  the  SNR’s  are  large  enough  such 
that  the  main  contribution  to  the  decorrelator  output  error  is  the  multiuser 
interference,  then  E{zkhk}  can  be  approximated  by, 

E{zkhk}  =  £'{2j;b*}(I- 2Pr(bi  in  error)) 

=  E{zkhk}{l-2Pek).  (19) 
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Equating  (19)  to  zero  and  using  (16)  together  with  (10),  we  get, 

E{zkhk]  =  E{pkkakkbk  +  pfAjbbfc  +  Ui  -  wjxjfc}bit  =  0,  (20) 

where  pkk  is  the  fcth  term  at  the  diagonal  of  P,  a^k  is  the  kih  diagonal  of  A, 
Pj.  is  the  ibth  column  of  P  with  its  fcth  element  deleted  and  so  is  Aj,. 

3.2.  Steady-State  Weights  and  Decorrelator  Output 

Solving  Eq.  (20),  we  show  in  the  Appendix  that  E{wjTc.kk>k}  =  A^pj,  We 
also  show  that,  Ejwf’xjtbfc}  =  AfcPfcWfc.  Therefore,  the  steady-state  weight 
resulting  from  the  decorrelating  algorithm  is  given  by. 

Wife  =  PJ  Vifc>  (21) 

Substituting  Eq.  (21)  in  (16)  with  (A-2),  we  drive  in  the  Appendix;  Eq.  (A-4), 
the  decorrelator’s  steady-state  outputs, 

Zk  =  ctkkiPkk  -  PkPk^rk)bk  +  nk-  pl’P^^nt.  (22) 

where  Pa:  is  obtained  from  P  by  deleting  the  ^th  row  and  column,  while  is 
the  feth  row  without  its  kth  element. 

IV.  PERFORMANCE  EVALUATION 


From  Eq.  (22)  for  user  one,  after  combining  the  noise  terms,  we  have, 

1  0 


zi  =  \/ar(l  -  pfPi^ri)6i -b  [  1  pJ 


[0  -Pr 


n. 


(23) 


Note  that  the  decorrelating  detector  perfectly  cancels  the  interfering  signal 
energy.  On  the  other  hand,  the  desired  user  bit  energy  is  modified  by  (1  — 
p^P^^n),  and  the  noise  variance  is  given  by 


where  p^  =  [  1  p^  ]  ^  ,  R  ^  = 

matrix  given  in  Eq.  (12). 

Let, 


—  Pt  ^  rt„X»'  Pc  n  ’ 

■  1  0 
.0  -Pr 

(l-pfP^^ri)2ai 


(24) 


and  R„  is  the  noise  covariance 


7i  = 


(25) 


’N 


The  BER  for  binary  PSK  as  a  function  of  SNR  is  given  by  (5(271),  where 
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V.  EXAMPLE  AND  RESULTS 


In  our  numerical  calculation  and  simulation,  we  assume  €2  =  0.4,  P21  =  0.2/y/0A,  pi2 
O.Q/y/0.6.  If  T2  <  T2,  then  Cj  <  ^2  (assumed  0.35)  and  hence  the  overlap  of 
$2  with  Si  is  larger  than  that  of  s^.  Therefore  we  take  p'12  =  0.25/\/0.35. 
With  the  same  argument,  p'21  =  0.55/-\/0.35r,pf'2^  =  \/ey^  -  ^0. 35/0. 4  and 

P22^  =  \/(l  -  f2)/(l  -  4)  =  \/0.6/0.65.  Finally,  pf^^  =  0.05/-y/(0. 65)0.4.  In 
summary. 


P^  = 


1 

0.3101 

0.9297 


0.3162 

0.9354 

0.0981 


0.7746 

0 

0.9608 


(26) 


Using  matrix  (26)  the  BER  performance  of  both  conventional  and  adaptive 
one-shot  detectors  is  depicted  in  Figs.  (3)  and  (4)  for  fixed  interference,  variable 
desired  user  energy  and  variable  interference,  fixed  desired  user  energy  respec¬ 
tively.  In  both  figures,  it  can  be  observed  that  especially  for  high  interfering 
energy  the  adaptive  decorrelator,  outperforms  the  conventional  one. 


Appendix 

It  is  easy  to  show  that,  =  0, 

E{iP^A.khk)hk}  =  AkPk,  E{nkhk}  =  0,  so  that 

E{iwjxk)hk]  =  AkPk-  (A-1) 

From  (10) 


x* 


[P^  r. 


■  A;t  O' 

b 

0  Okk 

.  . 

P'^Akhk  -t-  Tkakkhk  +  n*,. 


+  n* 


(A -2) 


Therefore, 


E{wlxkhk}  =  AkPk-Wk-  (A -3) 

From  Eq.  (16)  together  with  Eq.  (21), 

—  PkkO-kk^k  p'k  Akhk  Uk 

-PkPk'^iPl-^khk  +  rkOkkbk  +  nk)  (A  -  4) 

=  dkkiPkk  -  PkPk'^rk)bk +nk  -  plP^'^u.k. 
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SNR  1,  dB 


Fig.  4:  Performance  Comparison  of  .4daptivc  and  Conventional  DecoiTclator 
for  fixed  SNR  l  and  variable  Interference  valnes. 


Appendix  D 


Maximum  Signal-to-noise  Ratio  Data  Combining 
for  One-shot  Asynchronous  Multiuser  CDMA 

Detector  * 

Yeheskel  Bar-Ness  Nadir  Sezgin 


Abstract 

Using  a  one-shot  matched  filter  followed  by  a  decorrelating  detector, 
results  in  outputs  corresponding  to  the  main  user  and  the  left  (earlier 
bit)  and  right  (current  bit)  portion  of  the  code  of  the  other  users.  Due 
to  reduced  power,  these  outputs  will  result  in  a  degradated  probability 
of  error.  The  purpose  of  this  paper  is  to  propose  a  method  of  combining 
for  the  data  of  the  split  code  users.  We  suggest  to  use  combining  gain  in 
proportion  to  the  relative  user  delays.  The  resulting  decorrelating  detector 
always  performs  better  than  the  split-output  detector.  Such  diversity  data 
combining  eliminates  the  need  for  the  multiple  one-shot  filters,  one  for 
each  user. 


I.  INTRODUCTION 

CDMA  is  considered  to  be  a  promising  multiplexing  method  for  multiuser  per¬ 
sonal,  mobile  and  indoor  communications.  One  of  the  problems  a  designer  of 
such  a  system  is  faced  with  is  the  so  called  “near- far  problem,”  resulting  from 
excessive  Multiple  Access  Interference  (MAI)  energy  from  nearby  users,  com¬ 
pared  with  the  desired  user’s  signal  energy.  Power  control,  that  is,  adjustment 
of  transmitter  power,  depending  on  its  location  and  the  signal  energies  of  the 
other  users,  has  been  suggested  as  a  solution  to  this  problem.  But  it  requires 
a  significant  reduction  in  the  signal  energies  of  the  strong  users  i'n  order  for 
the  weaker  users  to  achieve  reliable  communication.  This  results  in  an  overall 
reduction  in  communication  ranges. 

An  optimum  near-far  resistant  multiuser  detector,  which  does  not  need 
power  control,  has  been  proposed  by  Verdu[l].  Its  complexity,  however,  is 
exponential  in  terms  of  the  number  of  users,  which  makes  it  unsuitable  for 
practical  situations.  Sub-optimum  decorrelating  detectors  which  are  based  on 
linear  transformation  of  the  sampled  match  filter  outputs  were  considered  in  [2] 

This  work  was  partially  supported  by  a  grant  from  Rome  Air  Force  Lab  (AFSC,  Griffss 
Air  Force  Base,  NY)  imder  the  contract  F30602-94-C-0118 
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for  the  synchronous  channel.  For  the  asynchronous  case,  Verdu  suggested  to  use 
a  one-shot  version  of  his  decorrelating  detector,  in  which  1  -|-  2(A'  -  1)  =  - 1 

{K  is  the  number  of  users)  filters  are  matched  to  one  full  code  of  a  specific  user 
(without  loss  of  generality,  called  the  first  user)  and  two  filters  matched  to  the 
each  other  user.  One  filter  matched  to  the  left  part  of  its  code,  corresponding 
to  times  between  zero  and  Tj  (the  delay  of  user  j  with  respect  to  the  first  user) 
and  to  the  right  part  of  the  code,  corresponding  to  times  between  Tj  and  T.  The 
decorrelating  detectors  were  also  proposed  as  a  first  stage  in  a  two-stage  mul¬ 
tiuser  CDMA  signal  separator  to  facilitate  tentative  decision  in  a  fixed  canceller 
stage[4-5]  assuming  knowledge  of  signal  energies,  and  in  an  adaptive  canceller 
stage  assuming  no  knowledge  of  energies[6-7].  One  may  easily  conclude  (see 
[7]  for  example)  that  with  the  one-shot  decorrelator/adaptive-canceller  we  can 
obtain  the  data  of  the  specific  user  to  which  the  one-shot  is  matched.  To  obtain 
the  other  users’  data  we  must  repeat  the  one-shot  matching  arrangement  to  the 
other  users.  Otherwise,  relying  on  the  left  or  right  portion  of  code  of  these  users 
will  result  in  reduced  performance. 

In  this  paper,  we  propose  proportional  ratio  diversity  data  combining  of  the 
left  and  right  portions  to  improve  performance,  instead  of  using  one  one-shot 
arrangement  for  each  user. 


II.  SYSTEM  MODEL 


For  the  asynchronous  multiuser,  the  equivalent  low-pass  signal  at  the  input  of 
the  matched  filter  bank  is  given  by. 


K 

^(0  =  +  iT  -  Tk)  +  n{t), 

k=zl  i 


(1) 


where  K  is  the  number  of  users,  ak,bk,  Sk  and  Tk  are  the  signal  amplitude,  user 
bit,  signature  waveform  and  the  relative  delay  of  ^th  user  respectively.  n(t)  is 
the  zero  mean  AWGN,  with  a  two-sided  power  spectral  density  of  No/2. 

For  the  sake  of  simplicity  we  will  restrict  ourselves  to  two  users  only.  Ex¬ 
tension  to  a  higher  number  is  relatively  simple. 

Representing  the  signal  of  Eq.  (1)  in  one-shot  of  the  ith  bit  of  user  one,  and 
without  loss  of  generality,  letting  i  =  0,  we  write, 

r{t)  =  +  y/a^^——s2{t)b2{—l) 

V^2 

+  \/a2(l  -  ez)  ..  ^  s§{t)b2{0)  +  n{t),  (2) 

y  L  —  €2  ■  ' 

where  0  <  r2  <  T  is  considered  unknown  or  estimated  with  an  error,  and 


4(0 


f  ®2(^  +  T  —  T2)  if  0  <t  <  T2 
1  0  if  ra  <  t  <  T 
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s^{i) 


£2 


Jo  if  0  <  t  <  r2 

\  S2(t  -T2)  if  T2<t  <T 


-\-T  —  T2^dt. 


(3) 


As  shown  in  Fig.  1,  we  apply  r{t)  to  a  bank  of  filters  by  using  as  a  second 
input  si{t),s^{t),  sf  (t),  respectively.  Note  that  we  have  assumed  the  knowledge 
of  r2.  The  result  of  this  paper  can  be  applied  to  the  case  where  r2  is  unknown 
or  has  been  estimated  with  error  (see  [8]).  Therefore,  without  loss  of  generality, 
T2  is  considered  to  be  known. 

The  output  of  the  first  filter. 


a:i(0)  —  y/aibi{0)  +  y/a^pi2f>2i—f) 

+  \/a2(l  —  e2)/!>f2^2(0)  +  ^i(O),  (4) 

I 

where  ni(0)  =  n(t)si{t)dt  is  a  zero-mean  Gaussian  variable  with  variance  of 
No/2,  and 

P12  =  s^it)si(t)dt 

P12  =  -7^=  I  (5) 

V  1  •"  ^2  Jo 


The  output  of  the  second  filter  is  given  by 

rT 


where 


*2(0) 


r(/)s2  {t)di 


^/d2  Jo 

■v/^^21^l(0)  +  Vd^^P22  b2{  —  f) 
+  y/  02(1  —  ^2)P22^b2{0)  +  1^2(0), 


^21  — 


P22 


P22  = 


and 

^2(0)  = 


T 

^lo  =  P12 

-  r{s^it)fdt=i 

C2  Jo 

1  /■’’ 


s^dt. 


(6) 


(7) 
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where  w*  is  the  ifeth  column  of  W  with  the  element  Wkk  deleted  and  Xk  is  the 
vector  obtained  from  x  by  deleting  Xk-  In  [9]  it  was  shown  that  such  W  can  be 
obtained  by  solving 

E{zkhk}  =  0,  k=  1,2,3-,  (15) 

where  hk  is  b  without  bk-  One  can  easily  show  that  if  the  signal-to-noise-plus- 
interference  ratio  at  any  Zj  is  large,  then  £'{zfcsgn(2:j)}  ps  E{zkbj}{l  —  Pej)  is  a 
high  SNR  approximation  of 

£;{2fcSgn(zjfc)}  =  0,  1,2,3.  (16) 

The  solution  of  Eq.  (14)  can  be  obtained  by  adaptively  controlling  the  weights 
Wmi  by  the  steepest  descent  algorithm,  which  simultaneously  reduces  the  ab¬ 
solute  value  of  the  correlation  between  the  outputs  of  the  decorrelator  and  the 
decision  on  all  other  outputs.  That  is,  the  weight,  wtk  is  controlled  by  the 
recursion, 

Wmi  <-  Wmi  -  n(zisgn{zm)),  l<l,m<3,  \i^m.  (17) 

In  [8]  it  is  shown  that  the  solution  of  Eqn.  (13)  is  given  by, 

Wfc  =  PJ  (1®) 

where  pj,  is  the  fcth  column  of  P  without  pkki  and  P*  is  the  matrix  P  without 
the  ibth  row  and  column.  Substituting  in  Eq.  (11)  we  get, 

zi(0)  =  -v/ar(l-pfPi  Vi)^i(0)  +  ”l>i(0)  I 

22(0)  =  •v/a2^(l  -  pI'PJ ^P2)^2(— 1)  +  ni>2(0)  (19) 

•2^3(0)  =  \/ 02(1  ~  f2)(l  “  pJPa  P3)^2(0)  +  ^£>3(0), 

where  np  =  [n£)i(0),  n£)2(0),  ud3(0)]^  =  Vn(0)  is  the  noise  vector  at  the  decor¬ 
relator  output. 

From  Eq.  (11)  we  notice  that  the  covariance  matrix  of  n£)(0)  is  given  by, 

Cn„  =  V^{n(0)n^(0)}V 

=  VPV^JVo/2.  (20) 

Hence  the  noise  variance  for  the  fcth  fictitious  user  (corresponding  to  the  left  or 
the  right  part  of  user  2)  at  the  decorrelator  output, 

<rDk  =  il-plP;^Pk)^  fc  =  l,2,3;  (21) 

where  we  used  the  fact  that  VP  is  a  diagonal  matrix  with  the  diagonal  entries 
(1  -  pJPj^  Vfc),  fc  =  1, 2, 3;  (see  Eq.(12).  Therefore, 

SNRfc  =  fc  =  1,  2,  3;  (22) 
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where  is  the  ifeth  diagonal  element  of  the  fictitious  users’  signal  amplitude 
matrix,  A. 


3.1.  Decorrelating  Output  Data  Combining 

If  we  use  ratio  combining  of  a  delayed  version  of  ^3(0)  with  ^2(0)  we  get  (see 

Fig.  2), 

22c(0)  =  22(0)  +  7%/!  -  e223(-l) 

=  ^/a^(l -P2P2  V2)fc2(-1) 

+  7A/a2(l  -  e2)(l  -  P3  P3  V3)^'2(-1) 

4- n/)2(0)  +  7«X)3(— 1).  (23) 

The  noises  71^32(0)  and  ni>3(— 1)  are  independent,  as  they  are  obtained  from  dif¬ 
ferent  time  periods.  Therefore,  by  using  Eqs.  (21)  and  (23)  we  get  the  combined 
output, 


SNRsc  = 


a2f2(l-P2P  V2)[H-Ty^ 


i-prp.7*p.o2 


fTj-prpj'ps 


[1  +  7' 


(24) 


One  can  show  that  7  =  y  result  in  maximum  SNR2C)  (see  Fig.  2). 

With  this  value  we  get  in  relation  to  uncombined  SNR2, 

SNR2c  =  (1  +  ^^^— 4?!^)SNR2.  (25) 

f2  1  P2  P2  P2 

In  relation  to  SNR3, 

SNR2C  =  (1  +  (26) 

1  £2  1  P3  P3  P3 

It  has  been  shown  in  [9]  that, 

0<l-pfPrVi=|^<l.  (27) 

where  |P|  is  the  determinant  of  the  cross-correlation  matrix,  and  P  and  |P,  |  is 
the  w'th  co-factor  of  P.  Since  0  <  e  <  1,  it  can  easily  be  concluded  that  SNR2C 
is  greater  than  both  uncombined  results,  for  any  e  and  any  correlation  matrix 
P. 

It  is  important  to  note  that  examining  equal  gain  data  combining  resulted 
in  an  improvement  only  for  certain  values  of  e. 


68 


IV.  RESULTS 

Using  the  linear  transformation  V  which  results  from  Eq.  (15),  then  applying 
(22)  for  the  uncombined  output  we  presented  in  Fig.  3,  the  probability  of  error 
as  a  function  of  SNR2  —  SNRj,  with  SNRi  =  8  dB  and  T2  =  0.4T.  The  Gold 
code  of  length  7  is  used  as  the  signature  sequences.  In  Fig.  4,  we  repeat  Fig.  3 
by  simulating  the  algorithm  of  (17)  instead  of  (15).  In  Fig.  5,  we  compare  the 
results  of  the  combined  data  for  the  second  user  obtained  with  (15),  versus 
simulation  with  (17).  We  notice  from  these  figures  that  the  probability  of  error 
of  the  combined  data  is  always  better.  Also  noticeable  is  the  fact  that  using  the 
bootstrap  algorithm  of  (17)  we  obtain  better  results,  particularly  at  low  SNR’s. 

Finally,  in  Fig.  6  we  compare  the  result  of  the  combined  data  (asynchronous) 
with  the  results  when  the  channel  is  synchronous. 
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Fig.  2:  Proposed  Combining  Scheme  for  Asynchronous  K  Users  at  bit  interval, 


Appendix  E 


Adaptive  Threshold  Setting  For  Multi-user 
CDMA  Signal  Separators  with  Soft  Tentative 

Decision  * 

Yeheskel  Bar-Ness  Nadir  Sezgin 


Abstract 

It  was  recently  shown  that  in  a  two  stage  decorrelator-canceUer  for 
multiuser  CDMA  separation,  soft  tentative  decision  can  help  to  improve 
performance,  particularly  when  interference  levels  are  moderate  to  low. 

The  threshold  points  of  the  soft  limiter  were  chosen  heuristically  and 
required  continuous  measurements  of  energies  of  all  signals  at  the  output 
of  the  decorrelator.  In  this  paper  we  propose  setting  these  thresholds 
adaptively  by  using  minimum  energy  at  the  outputs  of  the  canceller,  and 
hence  disposing  of  the  need  for  measurement.  The  resulting  decorrelator- 
canceUer  employs  two  algorithms;  one  for  the  canceUer  weights  and  the 
other  for  setting  the  threshold  points. 

I.  INTRODUCTION 

A  synchronous  multi-user  CDMA  receiver  that  employs  a  combination  of  a 
decorrelator  [1]  and  an  interference  canceller,  and  whose  weights  are  adaptively 
controlled,  was  presented  in  [2].  This  scheme  did  not  require  the  knowledge  of 
the  received  signal  energies,  but  performs  as  well  as  the  fixed  weight  counterpart, 
described  in  [3].  Both  used  a  hard  decision  on  the  decorrelator  output  to  obtain 
a  tentative  data  bit  estimate.  The  performance  of  this  two-stage  decorrelator- 
canceller  performs  very  well  when  the  interferer’s  signal-to-noise  ratios  are  high 
in  comparison  to  those  of  the  desired  signal.  The  performance  degrades  when 
the  former  becomes  the  same  level  or  lower  than  the  latter.  When  eliminating 
the  hard  limiter,  this  cascaded  two-stage  detector  performs  better  than  that 
with  the  limiter  at  a  low  interference-to-noise  ratio,  but  reaches  the  level  of 
the  decorrelator  performance  at  high  level  of  interference  (see  Fig.  1).  This 
motivates  the  use  of  a  soft  limiter  in  which  a  tentative  decision  is  made  only  at 
high  interference  energy.  Such  a  decorrelator-canceller  with  a  soft  limiter  was 
proposed  in  [4]  for  the  synchronous  case.  In  [5]  the  idea  was  extended  to  the 

This  work  was  partially  supported  by  a  griint  from  Rome  Air  Force  Lab  (AFSC,  Griffss 
Air  Force  Base,  NY)  under  the  contract  F30602-94-C-0118 
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asynchronous  case,  wherein  the  one-shot  decorrelator  of  Verdu[6]  was  used  as 
a  first  stage.  Zhang  and  Brady[7]  also  used  soft  decision  tentative  statistics  in 
their  multiuser  asynchronous  detector. 

The  question  that  remained  is  where  the  threshold  point  of  the  soft  limiter 
should  be  located.  In  all  the  aforementioned  references  the  threshold  value  was 
determined  heuristically  from  the  observed  values  of  the  decorrelator  outputs. 
If,  for  example. 


f  z/t  if  \z\  <  t 

\  sgn(2:)  otherwise  , 


(1) 


where  t  is  the  threshold  point  of  the  soft  tentative  decision  that  follows  the 


decorrelator.  In  [4],  t  was  chosen  such  that. 


tik  =  Plk 


[mk\}? 

E{Ui\} 


(2) 


where  Zk  is  the  output  of  the  decorrelator  corresponding  to  the  desired  user  and 
Z(  corresponds  to  the  interfering  /th  user,  pik  is  the  cross-correlation  between 
the  code  sequences  of  users  I  and  k.  A  similar  choice  was  used  in  [5],  for  the 
asynchronous  case.  Notice  that  for  setting  the  threshold  a  continual  measure  of 
energy  at  the  output  of  the  decorrelator  is  needed. 

In  this  paper  we  propose  an  adaptive  setting  of  this  threshold,  which  mini¬ 
mizes  the  detector  output  energy.  This  will  obviously  require  no  measurements. 


II.  SYSTEM  MODEL 


For  the  general  multi-user  the  equivalent  low-pass  signal  at  the  input  of  the 
matched  filter  is  given  by 


K  _ 

^(0  =  X!  X]  +  iT-Tk)  +  n{t),  (3) 

k  =  l  i 

where  K  is  the  number  of  users,  a*,  6*,  s*  and  Tk  are  the  signal  amplitude,  user 
bit,  signature  waveform  and  relative  delay  of  the  fcth  user.  n{t)  is  the  zero  mean 
AWGN,  with  a  variance,  cr^. 

For  the  synchronous  channel  encountered  in  down-link  communication  chan¬ 
nels  the  Tk  are  zero  for  all  k.  In  the  asynchronous  up-link  channels  the  rj,  are 
not  necessarily  equal  to  zero.  For  the  former  case  the  output  of  the  ifcth  matched 
filter  is  the  composite  of  bit  6*  and  all  interfering  bits  given  by  a  linear  com¬ 
bination  through  the  cross-correlation  factor  pkj  ■  In  matrix  notation  we  can 
write: 

a;(i)  =  PAb{i)  -f  n{i),  (4) 

where 

1  Pi2  •  •  •  PlK 
Pi2  1  •  ■  •  P2K 

(5) 

PlK  P2K  ■■■  1 
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and  A  =  diag(yd7,  •  •  • , -v/ok))  ^(0  =  [^*i(*)i '  ‘  ‘  i®  zero- 

mean  Gaussian  noise  vector  with  covariance  Rn  =  Po-.  For  the  asynchronous 
case  one  may  use  the  one-shot  matched  filter  proposed  by  Verdii  in  [6].  In 
this  case  Eq.  (3)  will  contain  a  linear  combination  of  hj{i)  and  hj{i  —  1) 
for  3  =  I...K,3  k.  That  is,  the  matrix  P  will  be  replaced  by  another 
square  matrix,  with  1  -|-  2(K  —  1)  =  2K  —  1  dimensions.  The  output  of  this 
one-shot  matched  filter  Xk{i)  will  contain  a  linear  combination  through  this 
(^2K  —  1)  X  (2K  —  1)  cross-correlation  matrix  of  bk{i)  and  bits  i  and  i  —  1  of 
all  other  users.  If  the  TkS  are  known  then  all  the  elements  of  this  matrix  are 
known  provided  one  knows  the  codes.  Therefore,  for  both  synchronous  and 
asynchronous  cases  the  output  of  the  matched  filters  is  given  by  Eq.  (4)  with 
the  corresponding  definition  of  b{i). 

For  simplicity  we  will  consider  the  two-user  case.  Extension  to  the  multiuser 
case  is  possible. 

Fig.  2  shows  a  block  diagram  of  the  proposed  two  stage  (decorrelator-canceller) 
detection  with  soft  decision  defined  in  Eq.  (1).  Let  A  =  1/t,  then 

(  1  if  A^  >  1 

f{z)  =  i  Az  if  -  1  <  Az  <  1  (6) 

(  —1  if  Az  <  —1. 

This  function  can  be  implemented  with  a  gain  whose  value  is  A  followed  by  a 
45°  soft  limiter  (see  Fig.  3). 

III.  DETECTOR’S  OPTIMAL  PARAMETERS  AND  ADAPTIVE 
CONTROL  ALGORITHM 

3.1  Optimal  Canceller  Weights 

Following  [1]  we  use  P~^  as  the  linear  transformation  for  the  decorrelator  to 
obtain 

z  =  Ab  +  i,  (7) 

where  ^  =  [^1,^2]^  is  a  Gaussian  vector  with  =  P~^a^ 

To  find  the  optimal  weight  for  the  canceller  stage,  we  write 


yi  =  Xi-  Wi2/(A2Z2). 


(8) 


Equating  the  derivative  of  E{yi}  with  respect  to  u;i2  to  zero  we  find. 


U'l2,opt(A2) 


E{xif{X2Z2)} 
E{P{X2Z2)}  • 


(9) 


Note  that  Eq.  (9)  describes  the  optimal  weight  in  terms  of  the  threshold,  I/A2, 
and  other  system  parameters.  Using  (A-1),  (A-3)  and  (A-6)  in  the  Appendix, 
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we  get 


U'l2,opt(-^2)  = 

P^/^{Gi{ai,  0:2)  +  ^2-^/^03(011,012)  +  Qjai)  —  Q(Q2)} 
G2(ao, «!,  «2)  +  +  «o)G3(q!i,  02)  +  Q(ai)  +  Q(a2)  ’ 


where 


Gi(Q'i,a2)  = 

02(010,  Oil,  a2)  = 

-  ^2) 

-  (2ao  +  ai)e“iG^ 

G3(q;i,Q!2)  = 

1-C?(ai)-Q(a2), 

ao  = 

Ol  = 

(1  - 

02  = 

(1  +  A2-ya^/(T,,2, 

TT  '^OC 

A1.WXV4.  XXi.XXAUN^X  jL>^.  Uj  ICttlll^  /\2 

infinity.  The  first  terms  in  the  numerator  and  denominator  clearly  go  to  zero. 
The  second  term  can  be  shown  to  go  to  zero  by  applying  L’  Hopital.  Thus  we 
are  left  with 

Sgjf-gfc; 

=  (11) 

the  same  result  obtained  in  [2]. 

3.2  Optimjil  Threshold  Level 
To  find  the  optimal  threshold,  by  using  Eq.  (8)  we  get 

r,  2  cifc/\  ^c(f('^2^2)^ 


dA2 


For  wi2  0 


(\  X  E{a;id/(A2Z2)/dA2} 

12(  2.opt)  ^{/(A2^2)d/(A2^2)/dA2}- 

Note  that  Eq.  (13)  describes  implicitly,  Aopt  as  a  function  of  the  weight  u;i2. 
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Using  Eqs.  (A-9)  and  (A-11)  of  the  Appendix  we  get 


W^12(A2,opt)  = 

p-,/a^[Gi  (01,02)  +  A2v^t?3(Qii.Q;2|)]  ^ 

^2(00,01, 02)  +  +  00)^3(01 , 02) 

Clearly,  for  this  definition  of  oo,oi  and  02  in  Eq.  (14),  we  mnst  use  A2,opt- 


3.3  Optimal  Threshold  Setting  for  Optimal  Detector 


If  we  define  the  terms  in  the  numerator  of  Eq.  (14)  as  A  and  those  in  the 
denominator  as  B,  then  using  U)i2,opt  and  A2,opt  in  both,  we  write  Eqs.  (10)  and 
(14)  as 


n^i2,opt(A2)  — 
n^i2(A2,opt)  = 


A  +  P^/a2[Q{o‘l)  -  Q(q^2)] 
B  +  Q(q;i)  4"  Q{o‘2) 

A 

B' 


(15) 


To  obtain  both  optimal  weight,  wi2,opti  and  A2,opt  we  must  equate  the  two 
equations  in  (15),  giving 


,,  ,  Pi/a2[Q(ai)  -  Q(«2)]  _  A 

tni2,opt(A2,opt)  -  g 

_  p^/^[Gi{ai ,  02)  +  A2-y/a2fj3(oi ,  ^2)] 
02(00,01,  02)  +  cr2^(l  +  0^)03(01, 02) 


(16) 


Solving  (16)  will  give  us  the  optimal  threshold  setting  A2,opt-  With  this  value. 


n^l2,opt(A2,opt) 


PV^[<?(ai,opt)  -  Q(a2,opt)] 

Q(oi,opt)  +  <3(02, opt) 


(17) 


where 


oi.opt  —  (1  ■”  A2,optV'n^/<^i)2i 

O2,opt  =  (1  +  A2,opf\/n^/(7,jj, 

and  (T,,2  =  A2(t/v/1  - 

From  Eq.  (15),  after  some  algebraic  manipulation  we  can  get  the  value  of 
the  optimal  threshold  ^2  =  I/A2: 

t2,opt  =  -7==(Q(ai)  -  <3(02))  (18) 

v^-p 

(1  +  5Ai?2,p)G3(oi,  02)  -  :^VSNR2,pGi{aua2) 

^104(01,02)  +  y/SNR2,p{Q{ai)  +  <3(o2))G3(oi,  02) 
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where  G4(ai,a2)  =  Q(ai)e-"'/2  -  Q(a2)e-“?/2  and  SNR2,p  =  02(1  -  p‘^)/a. 
This  equation  gives  implicitly  the  dependence  of  the  optimal  threshold,  topt,  as 
a  function  of  SNR2  and  p. 


3.4  Adaptive  Control  Algorithm 


The  following  two  steepest-descent  algorithms  simultaneously  control  the  weights 
Wij  and  the  gains  Xij,j  =  1,2, ..  .,K;j  ^  k  for  K  users  by  minimizing  the  en¬ 
ergy  E{y^i}.  That  is, 

E{yf(Xfj,wfj)}  (19) 

^  Va.,  E{yf{X’^j,wfj)}.  (20) 

IV.  PROBABILITY  OF  ERROR 


For  any  setting  of  A2,  the  optimal  weight,  tai2,opt(A2),  is  given  by  Eq.  (10).  With 
this  weight. 


))«(«■) 


“I  V  J-1  ^ 

/  ('72  —  A2i/ai^) 

•  ®^P( - 0-2 


V2 


+ 


V27r(7,,2  y_i  cr 


■)dv2 


72  +  y/dl 


)-Q( 


72  -  y/ai 


)] 


,  (^  +  A2-y/a2)^ , 

■  - 5  2 

^^’12 


where 


(21) 


=  P\/a2(c  -  1),  ^2  =  P^/a^ic  +  l), 

7i  =  Py/^{cr}2  -  1),  72  =  Py/d^{cT]2  +  1), 

and  c  =  Wi2{t2)/ Py/ai-  Note  that  for  <2  =  <2, opt,  using  Eq.'  (16)  we  get 

g  _  ^(Q'l.opt)  ~  Q(Q'2,opt) 

Q(«i.opt)  +  <3(02, opt) 


(22) 
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V.  NUMERICAL  RESULTS 


From  (19)  it  is  clear  that  the  optimum  threshold  setting  t2  is  independent  of 
SNRi.  The  solution  of  Eq.  (20)  numerically  as  a  function  of  p  is  given  in  Fig.  4, 
with  SNR2  as  a  parameter.  In  Fig.  5,  we  depict  <2, opt  as  a  function  of  SNR2  for 
p  =  0.7.  In  Fig.  6  we  compare  the  value  of  the  output  energy  obtained  with  the 
hard  limiter  and  with  a  threshold  calculated  according  to  the  heuristic  value 
of  Eq.  (1).  The  corresponding  probability  of  error  for  these  cases  are  given  in 
Fig.  7.  For  comparison  we  add  to  this  curve  the  probability  of  error  at  the 
output  of  the  decorrelator.  Fig.  8  depicts  the  same,  except  for  p  =  0.3. 

Note,  particularly  in  Fig.  7,  that  the  error  probability  with  the  heuristic 
setting  is  better  in  some  regions  than  with  the  optimum  setting  of  threshold. 
This  is  due  to  the  fact  that  the  optimization  was  performed  with  respect  to 
minimum  energy,  not  error  probability.  In  Fig.  8  the  performance  is  almost  the 
same.  Nevertheless,  the  advantage  of  adaptive  threshold  setting  is  in  disposing 
the  need  for  measurement. 

APPENDICES 

10 

£'{*1/(^2)}  =  E{{y/^bi  +  py/a2b2  +  ni)f(Tj2)} 

=  p^/a^E{b2f{V2)},  (■'^-1) 

where  we  used  the  fact  that  Z2  and  61  are  uncorrelated,  and  so  is  ni  with  Z2-  It 
is  easy  to  show  that  E{b2f{r]2)\b2  =  1}  =  E{b2f{'n2)\b2  =  —1}-  Therefore, 

E{b2f{ri2)}  = -)dm,  (A-2) 

y  J-00  >72 

where  with  62  =  li  >72  =  <^2^2  =  ^2(1/02  +  ^2)1  >72  =  ^2\/d^-  Separating  the 
integral  into  its  three  regions,  (-00,— 1),  [-1, 1]  and  (l,oo)  with  respect  to  772, 
then  changing  variable  x  =  (772  — 

1  r — 

E{62/(>72)}  =  ^[y_^ 

+  /  (cr,jX  +  A2V^e"®  f^dx] 

J  —  02 

=  -^(e-“^/2-e-“?/2)  +  A2V^ 

-v/27r 

•  [1  -  Q(ai)  -  Q(a2)]  +  Q(ai)  -  Q(«2), 

(A-3) 

where  ai  and  0:2  defined  in  Eq.  (10). 
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2.) 


E{f{m)}  =  :^[E{f\m)\h2  =  1}  +  E{f\r^2)\h2  =  -1}].  (A-4) 

It  is  easy  to  show  that  conditioned  on  62  =  1  and  62  —1  the  two  terms  are 

equal.  Therefore, 

E{f{ri2)}  =  E{f\rj2)\h^  =  1} 

-  -  /  /  (»?)exp( - ^ - )dr]2 

V27r<T„2  y_oo  2(7^^ 

1  rr.2.  ^  ,  (»?2-A2v^)\, 

+  /  )in,.  (A-5) 

Substituting  x  =  (t]2  —  A2v'a2)/o'»;2  in  the  first  integral  and  x  =  {t]2  +  \fE^I<^ri2 
in  the  second,  after  performing  the  integration  we  get 

E{/2(%)}  =  ^[(2ao-a2)e-“'/2 

—  (2q!o  +  «i)c  +  Oq) 

•  [1  -  Q(ai)  -  C?(a2)]  +  Qioci)  +  Q(a2),  (A-6) 

where  Oq,  ai  and  02  are  defined  in  Eq.  (10). 


=  E{i^b,  +  p^b2  +  } 

=  p^E{b2Z2^^},  (A-7) 

where  we  used  the  fact  that  Z2  is  uncorrelated  with  61  and  ni  and  that  f{i]2) 
is  a  deterministic  function.  It  is  easy  to  show  that  the  conditioned  on  both 
h  =  ±1,  -E'{&22^2d/(7;2)/d?72|62}  has  the  same  value.  Therefore, 

E{b2Z2^^  =  E{b2Z2^^\b2  =  1} 

ar)2  dr}2 


=  viTz  s  “p(-  )■'■> 


Changing  variable  x  —  {t]2  —  A2-y/ai')/a-,,2>  softer  integration,  and  using  (A-7),  we 
obtain 


d)^2  A2 

+  A2-v/a2[l  -  <5(02)  -  <3(a2)]}, 
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where  ai  and  are  defined  in  (10). 


S{/(A222)d/(A222)/dA2}  again  conditioned  on  62  =  ±1,  has  the  same  value. 
Therefore, 


=  E{f{r}2)z2 

aA2 

=  E{f{T]2)z2dfim)/dm\f>2  =  1} 


df iV2)- 
dr]2 


V^cr,,, 


/. 


m  I  iv2 
—  exp( - 


^2y/d^y 


)dr]2 


(A-10) 


Changing  variable  a;  =  (??2  -  A2v^)/<t,,  and  performing  the  integration  gives 


^{/(A2a2)-^3^}- 

^{-^[(2ao  -  02)6"“'/^  -  (2ao  +  ai)e"“'^^] 
A2  v27r 

+  <(1  +  «0)[1  -•«(«!)  -  QM)- 


5.) 

Pjj  =  i[Pr(ni  <  Wi2(t2)/(h2)  —  Py/d^b2  —  ^/oi) 

+  Pr{ni  >  W'i2(t2)/(h2)  —  P^/d^h  +  y/di)]- 

(A-12) 

By  using  the  fact  that  Pr{r}2  <  —1,62)  =  Pfijl2  >  Ij  ”^*2)  we  can  write, 

Pei  =  [^’K”i  <  '*^12(^2)  —  Py/a2^2  —  \/«T) 

+  Pr{ni  >  Wi2{t2)  —  P-v/^^2  +  \/^)]-PK’?2  >  1|^'2) 

+  'Wi2{t2)ll2  —  P^/E2^2  —  y/E\,  r]2  <  ll|) 

+  Pr(ui  >  ‘U>i2{t2)‘n2  —  P\fE2h2  +  -v/oi,  ??2  <  |1|)]- 

(A-13) 


By  letting 


=  Pyfadiic  -  1),  62  =  P'sfddiic  +  1), 

7i  =  Py/Eiic-f]  -  1),  72  =  P^/d^icn  +  1)> 
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Figure  1;  Performance  comparison  of  two-stage  detector,  with  Hard  Limiter  and 
without  Hard  Limiter;  SNRi  =  8  dB. 
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Figure  2:  Two-stage  (decorrelator/canceller)  Detector  with  Adaptive  Soft  Lim¬ 
iter. 


Figure  3:  Soft  Limiter  with  controlled  threshold. 


Figure  4:  Optimal  Threshold  as  a  function  of  correlation  coefBcient,  p,  for 
SNRi  =  8  dB. 
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Figure  5;  Optimal  Threshold  as  a  function  of  SNR2,  for  p  =  0.7. 


Figure  6;  Output  Energy  of  desired  user,  E{yl},  with  Hard  Limiter,  Heuristic 
Threshold  and  Optimal  Threshold;  SNRi  =  8  dB,  p  —  0.7. 
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Figure  7;  Bit  Error  Rate  performance  of  desired  user  for  Decorrelator,  Heuristic 
Threshold  and  Optimal  Threshold  Detector;  SNRi  =  8  dB,  p  =  0.7. 
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Figure  8:  Bit  Error  Rate  performance  of  desired  user  for  Decorrelator,  Heuristic 
Threshold  and  Optimal  Threshold  Detector;  SNRj  =  8  dB,  p  =  0.3. 
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Abstract 

A  stochastic  gradient-based  decorrelation  is  suggested  for  separation 
of  unknown  linear  mixture  of  signals.  It  is  shown  that  while  the  decorrela¬ 
tion  algorithm  is  similar  in  cost  to  the  LMS  algorithm,  its  rate  of  conver¬ 
gence  is  significantly  faster  making  it  more  attractive  for  signal  separation. 
Analysis  of  the  decorrelator  algorithm  shows  that  the  faster  speed  of  con¬ 
vergence  is  a  consequence  of  the  eigenvalues  spread  associated  with  the 
decorrelation  problem,  which  is  smaller  than  the  spread  associated  with 
the  corresponding  mean  square  problem.  Operation  of  the  algorithm  is 
illustrated  as  an  adaptive  multiuser  detector  in  a  multiple  carrier  CDMA 
system. 


1  Introduction 

In  this  paper  we  address  the  problem  of  the  recovery  of  unknown  independent 
sources  from  observations  of  a  linear  mixture  of  the  sources.  We  are  concerned 
with  a  multiple  input-multiple  output  system  where  each  output  is  an  unknown 
linear  combination  of  the  inputs.  The  following  linear  statistical  model  is  as¬ 
sumed: 

x=Ab-f-v  (1) 

where  x,  b,  and  v  are  random  vectors  with  values  in  Tl^ .  The  mixture  ma¬ 
trix  A  is  unknown  but  is  invertible.  The  problem  is  to  estimate  realizations  of 

*This  work  was  supported  in  part  by  a  grant  from  ROME  Air  Development  Center,  AFSC, 
Griffiss  Air  Force  Base,  NY,  under  contract  F30602-93-C-0118. 
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the  vector  b  given  observations  of  x.  This  type  of  problem  arises  in  numerous 
applications.  The  recovery  of  b  when  only  the  model’s  output  observations  x 
are  available  is  referred  to  as  blind  signal  separation.  To  find  a  solution  to  this 
problem  it  is  necessary,  in  addition  to  the  output  observations,  to  have  some 
information  about  the  statistical  properties  of  the  components  of  b.  Examples 
of  blind  signal  separation  applications  are  the  cancellation  of  cross-pol  inter¬ 
ference  in  dually  polarized  systems  and  multiuser  detection  in  multiple  access 
communications  systems.  The  latter  application  will  be  used  to  illustrate  the 
techniques  considered  in  this  paper. 

One  approach  to  signal  separation  is  to  view  it  as  an  interference  cancella¬ 
tion  problem,  where  at  each  channel  output  one  signal  is  considered  the  signal 
of  interest  and  the  other  signals  are  considered  interferences.  When  a  reference 
signal  is  available,  the  mean  square  error  (MSE)  between  the  output  and  the 
reference  signal  is  minimized  by  the  classical  Wiener  filter,  which  can  be  im¬ 
plemented  using  steepest  descent  algorithms  such  as  the  LMS  and  RLS.  Such 
an  example  can  be  found  in  [1],  where  the  LMS  algorithm,  operating  as  an 
interference  canceler,  is  used  as  a  separator  of  cross-polarized  signals. 

A  different  approach  is  to  treat  the  separation  problem  in  its  own  right. 
Notice  that  at  high  signal-to- noise  ratio  the  separation  problem  is  equivalent  to 
the  identification  of  the  system  A“^.  Beginning  with  the  early  eighties  a  class 
of  adaptive  signal  separators  have  been  proposed  that  in  effect  estimates  the 
inverse  of  the  mixture  matrix  A“^,  [2].  Their  research  was  mainly  aimed  at 
communications  applications,  such  ats  cross-pol  separation  in  dually  polarized 
links.  It  assumed  that  the  signals  to  be  separated  are  uncorrelated.  This  con¬ 
dition  however,  resulted  in  an  indeterminate  solution.  Separation  was  shown, 
provided  that  “some  knowledge”  is  available  by  which  the  signals  can  discrim¬ 
inated  (for  example  slightly  different  spectral  characteristics).  Later,  similar 
separator  structures  were  suggested  for  implementation  using  neural  networks 
[3].  They  solved  the  indetermination  problem  by  assuming  that  the  signals  to 
be  recovered  are  independent.  An  independence  test  was  approximated  by  way 
of  decorrelation  using  non-linear  functions  [4]. 

In  this  paper  we  analyze  the  performance  of  a  stochastic  gradient-based 
decorrelation  algorithm  for  blind  signal  separation.  The  algorithm  has  a  com¬ 
putational  complexity  comparable  with  the  LMS  algorithm,  but  is  shown  to  be 
faster  than  the  latter.  The  main  contribution  of  the  paper  is  to  prove,  through 
an  eigenvalue  analysis,  that  the  differences  in  the  speed  of  convergence  are  a 
consequence  of  the  different  eigenvalue  spreads  associated  with  the  implemen¬ 
tation  of  the  two  algorithms.  While  the  problem  statement  and  proofs  are  quite 
general,  the  specific  application  considered  is  a  synchronous  multiuser  multicar¬ 
rier  CDMA  communication  system  in  a  fading  environment.  In  this  application 
the  decorrelator  is  used  as  a  multiuser  detector.  Multiuser  detectors  are  a  sub¬ 
ject  of  continued  research  due  to  their  promise  of  increased  capacity  for  mobile 
communications,  personal  communications  and  other  wireless  networks.  The 
decorrelator  has  been  previously  suggested  as  a  low  complexity  near-far  resis- 
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tant  multiuser  detector  [5] .  Other  authors  proposed  non-adaptive  and  adaptive 
multiuser  detectors  incorporating  the  decorrelator  [6],  [7].  In  this  paper  we  focus 
on  the  transient  response  of  the  adaptive  decorrelator. 

The  system  model  is  presented  in  section  2.  In  section  3  we  define  separation 
criteria.  Adaptive  algorithms  and  the  convergence  analysis  of  the  decorrelator 
are  developed  in  section  4.  Section  5  contains  simulation  results.  Conclusions 
are  provided  in  section  6. 


2  System  Model 

The  general  system  model  is  given  in  eq.  (1)  where  the  received  signal  consists 
of  a  linear  transformation  of  the  transmitted  data  and  additive  noise.  To  be 
more  specific  and  to  lend  credibility  to  the  model,  we  consider  a  synchronous 
multiuser  multicarrier  CDMA  communication  system  such  as  depicted  in  [8]. 
The  synchronous  assumption  is  not  necessary  but  it  simplifies  the  notation.  An 
asynchronous  system  can  be  made  synchronous  by  oversampling  provided  all 
users’  signals  are  have  been  acquired.  The  model  presented  below  can  be  equally 
applicable  to  uplink  or  downlink  communications.  In  multicarrier  (MC)-CDMA 
an  information  symbol  is  transmitted  on  multiple  subcarriers.  The  subcarrier 
waveforms  are  chosen  mutually  orthogonal  (ex.  rectangular  RF-pulses  with 
frequencies  separated  by  an  integer  number  of  cycles).  Each  of  N  users  is 
assigned  a  unique  code  expressed  by  a  set  of  O/tt  offsets  of  the  subcarriers  phase. 
The  main  advantage  of  this  type  of  modulation  over  conventional  direct  sequence 
spread  spectrum  is  that  each  subcarrier  is  subject  to  frequency-nonselective 
fading,  thus  eliminating  the  need  for  the  complex  RAKE  receiver.  The  signal 
transmitted  by  the  n-th  user  for  the  i-th  bit  is  given  by; 

Sn(i)  =  V^bn{i)  Cn{m)p{t)  COS  (iOct  +  (2) 

m=:l 

where  is  the  waveform’s  energy.  The  information  symbols  are  assumed  binary 
£  {-1,1},  independent  and  equiprobable.  The  signature  code  sequence 
for  the  n-th  user  is  given  by  c„(m)  =  ±1/-\/M,  m  =  1, . . . ,  M,  where  M  is  the 
number  of  subcarriers.  p{t)  =  2/VT,  0  <  t  <  T  is  the  pulse  shape  and  T  is  the 
symbol  duration.  To  ensure  orthogonality  between  waveforms,  the  frequency 
offset  f2  is  chosen  to  be  an  integer  multiple  of  27r. 

The  functional  diagram  of  the  multicarrier-multiuser  CDMA  receiver  is  shown 
in  Figure  1.  The  n-th  user  signal,  after  coherent  amplitude  detection  with  each 
of  the  M  subcarrier  frequencies,  is  given  by  the  M-dimensional  vector: 

rn(f)  =  \/^^n(*)H„c„  (3) 

where  the  diagonal  matrix  H„  consists  of  the  channel  coefficients  hnm,'m  = 
1, . . . ,  M.  The  channel  parameters  h„m  assumed  to  be  Rayleigh  distributed. 
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The  paths  of  the  individual  subcarriers  are  assumed  to  be  independent.  Notice 
that  under  the  coherent  detection  assumption,  phase  changes  in  the  channel  do 
not  affect  the  output.  The  vector  cj  =  [cn(l))  •  •  - , c„(M)]  is  the  n-th  user  code 
sequence  vector.  Henceforth,  for  brevity  of  notation  explicit  time  dependency  is 
dropped  since  the  decision  on  (i)  requires  observation  of  signals  only  during 
that  interval.  Collecting  the  contributions  of  all  the  users  and  the  additive  noise, 
we  have  at  the  output  of  the  coherent  detectors: 

q  =  HEb  +  r]M  (4) 

where  E  is  a  diagonal  matrix,  E  =  diag(vO,  n  =  1, . . . ,  iV,  and  H  is  the 
matrix 

H  =  [HiCi,...,HjvCAf]  (5) 

b  is  the  vector  of  information  bits,  and  the  components  of  the  noise  vector 
riM  are  assumed  i.i.d.  zero-mean  Gaussian  with  variance  cr^.  The  signal  after 
cross-correlation  with  the  code  of  the  j-th  user  can  be  written, 

xj  =  cjHEh  +  cj  TIM  (6) 

Finally,  collecting  the  terms  for  all  the  users  the  received  signal  can  be  expressed 
in  vector  form, 

X  =  AEb  -b  V  (7) 

where  A  =  C^H  is  the  mixture  matrix,  and  v(i)  is  the  noise  vector  with  the 
j-th  component  equal  to  cj t}m-  Since  the  channel  coefficients  are  unknown,  the 
matrix  A  is  unknown.  It  can  be  readily  shown  that  the  noise  covariance  matrix 
is  given  by, 

E  [vv^]  =  (8) 

Notice  that  the  MC-CDMA  model  in  eq.  (7)  is  the  same  as  the  general  model 
assumed  in  (1).  Multiuser  detection  has  commonly  been  considered  for  uplink 
communications.  With  direct  sequence  CDMA,  orthogonality  between  the  sig¬ 
nature  waveforms  is  maintained  in  the  synchronous  uplink  channel.  However, 
the  same  does  not  hold  for  MC-CDMA,  where  subcarriers  fade  independently. 
Hence,  decorrelation  cannot  be  achieved  by  a  simple  inversion  of  the  signature 
waveforms  correlation  matrix.  Analysis  of  the  convergence  rate  of  adaptive 
decorrelation  in  the  uplink  channel  is  made  very  difficult  by  the  fact  that  the 
matrix  AE  is  not  symmetric.  Numerical  examples  in  section  5  consider  the 
uplink,  however  the  analysis  in  section  4  assumes  that  AE  is  symmetric.  This 
model  represents  the  downlink  of  a  MC-CDMA  system.  In  this  case,  multiuser 
techniques  could  be  useful  to  improve  performance  at  the  mobile.  Systems  re¬ 
quiring  multiuser  detection  are  expected  to  operate  near  capacity  most  of  the 
time,  thus  the  mobile  could  operate  assuming  all  users  are  present.  All  signals 
transmitted  by  the  base  station  and  received  at  the  mobile  experience  the  same 
channel,  i.e.,  H„  =  H.  With  this  condition,  the  mixture  matrix  is  given  by 
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A  =  C^HC,  and  is  clearly  symmetric.  Furthermore,  it  can  be  assumed  that  all 
signals  are  received  at  the  mobile  with  equal  energies  ,  i.e.,  E  =  I,  where  I  is 
the  identity  matrix  (a  constant  gain  common  to  all  users  can  be  incorporated 
in  the  matrix  A).  The  model  for  the  received  signal  becomes, 

X  =  Ab  +  V 

where  the  mixture  matrix  A  is  symmetric  and  unknown.  Direct  detection  of  the 
elements  of  vector  x  is  impeded  by  the  presence  of  the  cochannel  interference. 
Therefore,  the  goal  is  to  add  a  receiver  stage  that  provides  signal  separation 
through  cochannel  interference  cancellation. 


3  Signal  Separation  Criteria 

Two  signal  separation  criteria  are  considered:  separation  through  minimization 
of  the  Mean  Square  Error  (MSE)  and  separation  through  signal  decorrelation. 
The  criteria  are  used  to  develop  control  algorithms  for  the  network  weights.  The 
network  weights  are  represented  by  a  matrix  W  =  [wi , . . . ,  w^r] ,  where  Wn  is  the 
N  xl  weight  vector  applied  to  the  network  input  vector  x  to  obtain  the  output 
for  the  n-th  user,  j/n  =  w^x.  An  estimate  of  the  data  bits  is  then  obtained  from 
b  =  sgny,  where  y  =  W^x.  The  following  conditions  are  assumed  to  hold  for 
the  data  and  the  channel:  (1)  the  components  of  the  data  vector  b  are  assumed 
independent,  binary,  and  equiprobable,  (2)  the  mixture  matrix  is  symmetric,  (3) 
the  mixture  matrix  A  is  strictly  diagonally  dominant,  i.e.,  an  >  Y^j^i  |aij|  for 
i  =  1, . .  .,N.  The  latter  condition  ensures  that  in  the  absence  of  noise,  a  correct 
decision  is  made  with  respect  to  each  of  the  desired  users  with  probability  one. 
This  condition  also  implies  that  the  matrix  A  is  invertible  [9,  p.349]. 


3.1  MSE  Separation 

The  MSE  signal  separator  minimizes  the  mean  squared  error  between  its  output 
and  a  reference  signal.  Typically  the  reference  is  initially  supplied  by  a  training 
signal.  When  the  adaptive  weights  converged  and  the  errors  with  respect  to  the 
training  signal  are  small,  the  detector  is  switched  to  operate  in  decision  directed 
mode,  and  the  reference  signal  is  supplied  by  the  estimated  symbol.  The  MSE 
signal  separator  is  the  optimal  linear  receiver.  The  network  weights  are  given 

by, 


=  arg  min  E 


^  Rr^Ri 


b- W^x 


(9) 


bx 


where  R^,  =  E  [xx^]  is  the  correlation  matrix  of  the  inputs,  and  Rj^  =  E  j^xb^ 
is  the  cross-correlation  matrix  between  the  estimated  bits  and  the  input.  We 
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will  show  that  the  application  of  this  weights  matrix  leads  to  signal  separation 
in  the  high  SNR  case.  Neglecting  the  noise  contribution  we  have: 


~  Ai;[bb^]A^ 

=  A2 


(10) 


where  we  used  E  =  I,  since  the  signals  are  assumed  independent  between 

users.  The  cross-correlation  matrix  is: 


=  s[(Ab  +  v)b’'] 

=  AP 


(11) 


where  P  is  a  diagonal  matrix,  {P};^  =  1  —  2Pr  p,-  9^  fcj 
shown  in  the  appendix  that: 


E  bibj 


1  2Pet  *  —  j 
0  i  ^  j 


1  —  2Pei,  and  it  is 
(12) 


The  weights  are  given  by  =  A“^P.  When  the  effect  of  the  noise  can  be 

neglected,  this  linear  transformation  applied  to  the  input  vector,  recovers  the 
transmitted  signals  vis.,  y  =  =  Pb  +  ~  Pb. 


3.2  Signal  Separation  by  Decorrelation 

The  decorrelator  seeks  to  decorrelate  each  output  y„  from  all  the  other  users 
bit  estimate.  This  criterion  can  be  expressed  as: 

Wd  =  arg  |by^  =  A  j  (13) 

where  A  is  a  diagonal  matrix.  We  now  develop  an  expression  for  the  weight 
matrix  solution  to  the  decorrelation  criterion  of  eq.  (13)  and  then  show  that 
this  weight  vector  indeed  performs  signal  separation.  Again  assuming  negligible 
noise  and  using  the  results  of  the  appendix,  we  have: 

E  by^  ~  E  Jbb^AW  ,  . 

and  from  eq.  (13)  the  decorrelation  criterion  is  equivalent  to  the  condition 
PAW  =  A.  It  results  that  the  decorrelation  weight  matrix  solution  is  given  by: 

Wd=A-^Ai  (15) 

where  Ai  =  AP  is  a  diagonal  matrix.  Indeed  when  this  weight  matrix  is 
applied  to  the  noiseless  input  vector  x  it  results  in  signal  separation  vis.,  y  = 
Wjx  ~  WjAb  =  Alb.  Note  that  while  the  amplitude  of  the  separated  signals 
cannot  be  determined  by  this  procedure,  the  bit  detection  is  not  affected  by  this 
indetermination  since  signAib  =signb. 
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4  Adaptive  Decorrelation 

Adaptive  signal  separation  algorithms  can  be  directly  derived  from  the  MSE 
and  decorrelation  criteria  formulated  in  eqs.  (9)  and  (13),  respectively.  The 
MSE  criterion  can  be  implemented  using  the  well  known  LMS  adaptation.  In 
this  section  we  formulate,  analyze  and  compare  the  decorrelation  algorithm  with 
the  LMS  algorithm.  LMS  update  of  the  signal  separation  network  weights  is 
expressed  by  the  relations: 

W(A;+1)  =  W  (k)  +  fix{k)e^{k)  (16) 

e{k)  =  b(ik)-W^(fc)x(fc) 

b  (k)  =  signb  (k) 

The  convergence  analysis  of  the  LMS  algorithm  is  well  known.  Under  a  set  of 
assumptions  often  referred  to  as  independence  assumptions  ,  a  necessary  and 
sufficient  condition  for  convergence  is  that  the  step  size  parameter  meets  the 
condition  0  <  fJ.  <  2/Amax(Rr) ,  where  Amax(R.i)  is  the  largest  eigenvalue  of 
the  covariance  matrix  R® .  It  is  also  well  known  that  the  convergence  speed  of  the 
LMS  algorithm  is  determined  by  eigenvalue  spread  x  =  Amax(Ra:)  /Amin  (Rx)  • 
A  stochastic  gradient-based  decorrelation  algorithm  can  be  formulated  di¬ 
rectly  from  eq.  (13).  It  seeks  to  null  the  instantaneous  cross-correlation  between 
the  outputs  of  the  signal  separation  networks  and  the  detected  symbols: 

W{k  +  l)  =  W{k)-ti  h{k)y'^{k) -dmg(h{k)y'^{k)^  (17) 

diag(W(A  +  l))  =  I 

b  (k)  =  signb  {k) 

where  the  operator  diag  generates  a  diagonal  matrix  of  the  elements  on  the 
main  diagonal  of  the  argument,  and  W  (k)  =  W(I:)-diag(W(A;)).  Note  that 
the  adaptive  decorrelator  has  a  computational  complexity  similar  to  the  LMS. 
Signal  separation  by  decorrelation  is  different  from  the  MSE  separator.  The 
MSE  weight  vector  for  any  user  is  derived  using  previous  symbol  estimates  of 
that  user  as  a  reference  signal.  In  contrast,  the  decorrelator  can  be  interpreted 
as  using  reference  signals  derived  from  all  other  users  bit  estimates.  Suggesting 
an  operation  where  processing  for  each  user  is  helped  by  the  other  users,  this 
structure  has  been  referred  to  as  bootstrapped,  [10] 

4.1  Analysis  of  the  Adaptive  Decorrelator 

The  LMS  decorrelator  have  been  observed  to  converge  faster  than  the  conven¬ 
tional  LMS  error  algorithm  [10].  We  prove  that  the  eigenvalue  spread  for  the 
adaptive  decorrelator  is  smaller  than  for  the  LMS,  thereby  providing  the  rea¬ 
soning  for  the  faster  convergence  speed  of  the  former. 
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Define  the  following  matrix  and  vector  partitions:  wf  =  ^ 


adaptive  decorrelator  in  eq.  (18)  can  be  rewritten  as  follows: 


[w  (A:  +  1)]  =  [I  -  fiRg-]  E  [w  (fc)]  -  [B  {k)  x,  {k)  w,  (k)]  (21) 

where  =  E  [Bx^] .  Form  this  relation  we  deduce  that  the  convergence 

properties  of  the  adaptive  decorrelator  are  controlled  by  the  eigenvalues  of  R;g^- 

Proposition  1.  For  convergence  of  the  decorrelation  algorithm  the  following 
condition  has  to  be  satisfied: 

0  A*  ■'i  2  /  Ajjiax  (Rfl®)  (22) 

Proof  The  proof  follows  from  the  well  known  convergence  condition  of  the 
LMS  algorithm  [11],  and  the  analogy  of  the  LMS  and  decorrelator  algorithms. 
Proposition  2.  The  following  relation  exists  between  the  largest  eigenvalues  of 
and  Rj-^-’ 

Amax  (Rjj;)  >  -^max  (23) 
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Proof.  The  matrix  R;g^^  can  be  expressed  as  a  partition  of  the  matrix  Rj-^; 


ri 


R-nx  / 


(24) 


Define  the  ratio 


Piz) 


z^R-  z 

OT 


(25) 


Applying  the  Courant-Fisher  Minimax  theorem  [12],  A^ax  (%j:)  =  Pi^)  ■ 
Using  the  matrix  partition  in  eq.  24  and  maximizing  p{z)  over  the  restriction 
zi  =  0,  where  zi  is  the  first  component  of  the  vector  z,  we  get 


Amax  (R^,,)  >  maxp(z)  =  Aroax  (Rfij) 


(26) 


•  1=0 


Proposition  3. 
Rjg-  and  R-j-^; 


The  following  relation  exists  between  the  smallest  eigenvalues  of 


Amin  (^Bx) 


(27) 


Proof.  According  to  the  Interlacing  property  which  follows  from  the  Courant- 
Fisher  Minimax  theorem  [12],  and  noting  that  R^j  is  the  {N  —  1)  x  {N  —  1) 
leading  principal  submatrix  of  Rj^,  we  have: 


Amin  —  ^N-1  (l^Bx)  “  (^Bx) 

Proposition  4.  If  the  noise  is  negligible  then  the  eigenvalues  of  Rj,  are  equal  to 
the  square  of  the  eigenvalues  o/R-g-^. 

Proof  We  have  the  relation  R-g^  =  PA  ~  A,  and  for  negligible  noise  we 
have  Rj;  ~  AEEA^.  Consequently,  Rj,  ~  RJ^R-j-^.-  Since  R-^^  (with  the 

approximation  P  ~  I,  is  symmetric,  it  follows  that  Aj  (Rj,)  =  A?  for 

1  <  i  <  A.  In  particular,  we  have  for  the  eigenvalue  spread. 


x(R-x) 


■^max  (Rx) 
^min  (Rx) 


(29) 
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Proposition  5.  The  eigenvalue  spread  associated  with  the  LMS  algo¬ 
rithm,  is  larger  than  the  spread  x  associated  with  the  decorrelator  algo¬ 

rithm. 

Proof.  From  Propositions  2,3  and  4  we  have: 


=  (30) 

Proposition  5  provides  the  explanation  why  the  decorrelator  algorithm  is 
faster  than  the  LMS  algorithm. 

Proposition  6.  The  following  relation  exists  between  the  largest  step  size  pa¬ 
rameters  of  the  decorrelator  and  LMS  algorithms: 


The  LMS  and  the  decorrelator  algorithms  were  compared  in  terms  of  the  prob¬ 
ability  of  error  for  detecting  the  transmitted  data.  Consider  the  probability  of 
bit  error  for  the  i-th  user.  Since  the  data  bit  can  take  on  one  of  two  equiprob- 
able  values  the  total  probability  of  error  is  equal  to  the  conditional  probability 
P  (bi  ^  bi  I  bij  .  The  probability  of  error  at  the  output  of  the  separator  can  be 
computed  as  a  function  of  the  weight  vector  wf : 

Pei  =  P  (bi  ^  bi  I  bi^  (33) 

=  P  (wfx  >  0  I  6,- = -l) 

=  [-P(wfx  >  0  I  = -l)] 

6<=-l 

=  2^-^  E  ^ 

bi=-l 

2 

where  Q{z)  =  e~  V  dy,  t^j  is  the  j-th  element  of  the  vector  tj  =  wf  AE. 

The  convergence  of  the  algorithms  was  studied  by  simulations.  The  simula¬ 
tions  consisted  of  a  iV  =  4  users  system  each  using  an  M  =  15  chip  Gold  code 
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spreading  sequence.  The  multipath  fading  channel  was  modeled  as  Rayleigh 
random  variables  independent  between  carriers.  The  fading  characteristics  were 
assumed  fixed  through  each  simulation  run  and  independent  between  runs.  The 
probability  of  error  was  evaluated  with  respect  to  the  first  user.  The  SNR  is 
defiend  with  respect  to  the  first  user.  The  signal-to-interference  ratio  (SIR)  is 
defined  as  the  power  ratio  of  the  designated  user  to  each  of  the  other  users  (as¬ 
sumed  to  have  equal  power).  Figure  2  shows  the  average  learning  curves  of  the 
probability  of  error  for  SNR  =  12  dB  and  SIR  =  -3  dB.  Shown  is  the  average 
of  20  runs  (20  independent  channels).  The  channel  was  held  fixed  during  each 
run.  This  simulation  represents  the  uplink  channel  model.  Signals  are  received 
through  independent  channels,  and  the  mixture  matrix  is  asymmetric  and  given 
by  A  =  C^H,  where  H  is  given  by  eq.  (5).  The  faster  convergence  speed  of  the 
decorrelator  is  strikingly  evident.  To  evaluate  the  steady  state  response  of  the 
algorithms,  the  following  error  probability  benchmarks  are  shown  in  Figure  3; 
(1)  decorrelation  with  (C^C)”^  (channel  effects  ignored),  (2)  decorrelation  with 
A“^  (known  channel),  and  (3)  single  user  (no  multiuser  inferference).  The  sce¬ 
nario  was  the  same  as  in  the  previous  figure,  except  that  Figure  3  represents  a 
single  run  generated  with  the  same  channel  as  the  one  used  to  compute  bench¬ 
mark  (2).  The  mixture  matrix  A  for  this  run  is  given  by; 


.0.8605 
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Notice  that  decorrelation  using  the  inverse  of  the  signature  waveforms  correla¬ 
tion  matrix  (curve  9^11),  results  in  considerable  performance  degradation.  The 
simulation  results  clearly  show  that  the  adaptive  decorrelator  has  a  higher  rate 
of  convergence  than  the  LMS  algorithm,  as  well  as  converging  closer  to  the  error 
for  known  channel  (curve  #2).  Figure  4  shows  the  average  of  20  runs  for  the 
case  of  symmetric  mixture  matrix.  This  is  the  downlink  case  with  SIR  =  0  dB 
and  all  users  experiencing  the  same  channel.  The  decorrelator  has  faster  conver¬ 
gence  rate  than  the  LMS,  but  the  difference  is  smaller  than  in  the  asymmetric 
case.  This  is  due  to  the  fact  that  negative  SIR  actually  helps  the  decorrelator 
which  is  using  decisions  of  the  other  users. 

To  estimate  the  convergence  region  of  each  algorithm  we  took  the  approach 
suggested  in  [13].  A  figure  of  merit  7  is  defined  that  relates  the  initial  and  final 
probabilities  of  error,  Pei,  and  Pej,  follows: 


7  =  1-^  (34) 

when  Pef  <  Pein,  7  RS  1.  Note  that  7  =  0  corresponds  to  no  convergence  is 
indicated  by  Pef  -  Pein-  Figure  5  shows  the  convergence  curves  when  SIR  = 
-3  dB,  and  the  SNR  is  varied  from  -6  dB  to  -|-14  dB,  to  result  in  the  initial 
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probabilities  of  error  indicated  on  the  abscissa.  The  figure  shows  that  not  only 
is  the  decorrelator  faster,  but  it  also  has  a  wider  region  of  convergence  than  the 
LMS. 

6  Conclusions 

In  this  paper  we  analyzed  the  transient  response  of  a  stochastic  gradient-based 
decorrelation  algorithm  and  showed  that  it  has  a  faster  convergence  rate  than 
due  to  a  lesser  eigenvalue  spread.  Furthermore,  the  algorithm  was  shown  to 
have  wider  regions  of  convergence  than  the  LMS.  The  decorrelation  and  the 
LMS  algorithms  are  of  the  same  complexity.  The  application  of  the  adaptive 
decorrelator  was  suggested  as  a  multiuser  detector  for  a  synchronous  multiple 
carrier  CDMA  system  in  a  fading  environment. 

A  Evaluation  of  E  [6*6*] 

This  appendix  contains  the  evaluation  of  E  6, A,]  .  Prom  the  definition  we  have, 

■S' =  Y}^.  iM^iP  (%^i) 

\  .  I  (35) 

The  sum  over  6,-  can  be  separated  into  correct  and  erroneous  decisions, 

=  0.5x2(p(?,  =  6,|6,)+p(T,^6,|t,))  (36) 

=  l-2Pei 

where  P(e,)  is  the  probability  of  error  for  the  i-th  user. 


(37) 


(38) 


B  Evaluation  of  E  [6*6j] 

From  the  definition  we  have, 

E  [6i6j]  =  Eb-^h-bibiP  (bibj^ 

The  joint  probability  of  bi  and  bj  can  be  developed  as  follows: 
P{bibj)  =  P(6.|6i)i^(6*) 

=  J:i,P(bi\bi,bj)p{bi)P{bj) 

=  0.25(^p(bi  =  bi\bj)+p(bi^bi\bj)) 
=  0.25(l-2Pei) 
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Thus  we  get, 

=  0.25(l-2Pei)Ej.,i,'^i^i 

=  0 

The  results  in  eqs.  (36)  and  (39)  are  combined  and  shown  in  eq.  (12). 


E 


bi  bj 


(39) 
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Figure  1:  Multicarrier-multiuser  CDMA  receiver. 


Figure  2;  Learning  curve  of  the  average  probability  of  error  of  the  first  user. 
Asymmetric  mixture  matrix.  SNR  =  12dB,  SIR  =  -3dB. 
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Figure  3:  Learning  curve  of  the  average  probability  of  error  for  a  single 
Asymmetric  mixture  matrix.  SNR  =  12dB,  SIR  =  -3dB. 


run. 
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Figure  4:  Learning  curve  of  the  average  probability  of  error.  Symmetric  mixture 
matrix.  SNR  =  12dB,  SIR  =  0  dB. 
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Figure  5:  Convergence  regions  after  100  iterations  for  the  probability  of  error 
of  the  first  user.  Asymmetric  mixture  matrix  SNR  =  12dB,  SIR  =  -3  dB. 
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Abstract 

In  this  work,  we  present  a  new,  computationally  simple  scheme  (termed  the  multi-shot  approach)  to 
separate  and  detect  multiuser  signals  in  asynchronous  CDMA  communication  systems.  By  exploring  the 
structure  inherent  in  the  matrix  decomposition  of  properly  arranged  data,  obtained  from  multi-shot  matched 
filtering,  we  propose  a  near-far  resistant  multiuser  detector.  The  proposed  multiuser  detector  combines  the 
multi-shot  matched  filter  outputs  through  a  matrix  filtering  and  de-biasing  processing,  suggested  by  the 
inherent  structure  in  data.  Simulation  results  support  the  newly  proposed  detection  scheme. 


I  Introduction 

Code  Division  Multiple  Access  (CDMA)  is  an  accepted  schenae  for  future  high  capacity  digital  wireless/ cellular 
communications.  The  need  for  a  large  number  of  users  within  a  finited  bandwidth  dictates  the  implementation  of 
a  semi-orthogonal  set  of  signature  waveforms,  such  as  Gold  or  Kasami  codes.  The  cross  correlation  between  the 
user’s  code  will  cause  Multiuser  Interference  (MI).  A  conventional  digital  receiver  that  ignores  the  existence  of 
other  users  will  degrade  quickly  with  the  increase  of  interference.  Further,  the  existence  of  many  users  will  heavily 
interfere  with  the  desired  weak  user.  This  is  referred  to  cis  the  ‘near-far  problem’  in  CDMA  systems.  Hence  a 
multiuser  detector  is  necessary  for  multiple  access  (MA)  communication  systems.  Many  near-far  resistant  mul¬ 
tiuser  detectors  have  been  proposed  and  analyzed  for  synchronous  CDMA  systems  [1,2,3].  For  the  asynchronous 

tThis  work  was  supported  in  part  by  Rome  Air  Force  Lab  under  contract  F30602-94-C-0135,  and  by  the  Office  of  Sponsored 
Research,  NJIT. 
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case,  that  is,  when  the  received  users’  signals  are  not  synchronized  in  time  such  as  in  uplink  (from  mobiles  to 
base  station)  wireless  communications  [4,5],  the  idea  of  a  ‘one-shot  decorrelator’  was  proposed,  where  a  bank 
of  filters  is  matched  to  one  of  the  users’  code,  while  the  others  are  matched  to  the  left  and  right  parts  of  other 
users  code.  This  merely  increases  the  dimension  of  the  matched  filtering  output  vector  and  the  cross-correlation 
matrix.  Except  for  the  large  dimension,  the  near-far  resistant  multiuser  detector  for  an  asynchronous  channel  is 
similar  to  the  synchronous  case.  It  was  recently  demonstrated,  however,  that  for  some  relative  delays  between 
the  user  codes,  the  correlation  matrix  becomes  singular,  so  that  the  above-mentioned  linear  transformation-based 
method  failed. 

Instead,  in  this  work,  we  present  a  new  multi-shot  approach  to  separate  and  detect  multiuser  signals  in 
asynchronous  CDMA  communication  systems.  The  proposed  method  first  explores  the  structure  inherent  in  the 
matrix  decomposition  of  the  properly  arranged  data,  obtained  from  multi-shot  matched  filtering.  Then  we  make 
a  matrix  approximation  to  make  the  tentative  multiuser  information  bit  estimation  trackable  and  simple.  After 
analyzing  the  bias  intorduced  by  matrix  approximation,  we  further  proposed  a  matrix  filtering  and  debiasing 
scheme  to  improve  the  overall  detection  performance.  Simulation  results  support  the  newly  proposed  detection 
scheme. 

II  Problem  Formulation 

In  CDMA  communication  systems,  the  data  r(t)  at  the  receiver  is  actually  a  mixture  of  K  multiuser  signals,  that 
is, 

K 

fc=i  » 

where  K  is  the  total  number  of  users,  r*,,  Sk{t),  and  hk{i)  are  the  amplitude,  transmission  delay,  signature 
waveform  (of  duration  T),  and  information  bit  (at  the  ith  symbol  interval)  of  the  /feth  user,  respectively.  'w{t)  is 
an  additive  white  Gaussian  noise  process  with  zero  mean  and  a  two-sided  power  spectral  density  (PSD)  of  cr^. 

For  the  synchronous  transmission  channel,  the  transmission  delays  associated  with  all  the  users,  r^’s,  are  equal 
and  can  be  treated  as  zeros.  We  showed  in  [6]  that  the  outputs  of  a  bank  of  K  matched  filters,  sampled  at  the 
end  of  the  zth  sybmol  interval  {t  =  iT)  form  the  sufficient  statistic  of  all  the  users’  information  bit,  {6fc(i)}^^i, 
within  the  zth  symbol  interval.  Hence,  the  one-shot  approach  is  an  appropriate  way  to  detect  and  separate 
multiuser  information  bits  [1,2,3].  However,  for  an  asynchronous  transmission  channel,  situation  becomes  more 
complicated,  since  delays  of  all  the  users,  r^’s,  are  not  the  same  in  general.  In  the  following  sections,  we  propose 
a  multi-shot  approach  and  associated  signal  processing  techniques  to  get  a  near-far  resistant  multiuser  detector 
through  matrix  decomposition,  approximation,  and  debisaing  processing. 

Without  loss  of  generality,  we  make  the  assumption  that  delays  of  all  the  K  users  satisfy  the  relationship 
0  <  T"!  <  72  <  •  •  •  <  rK  <  T.  Recall  that  the  signature  waveforms  of  all  the  users,  Sfc(t)’s,  are  of  duration  T, 
therefore,  any  delay  values  larger  than  the  symbol  interval  T  can  be  treated  into  another  bit.  We  then  consider 
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the  output  of  the  fcth  matched  filter,  sampled  at  the  t  =  iT  (i  =  0,  1,  . . .)  instant, 

f(!  +  l)T+T)i 


xk{i}  =  / 

Ji  T  +  TJ, 


r{t)  Sk{t  -  iT  -  Tk)d.t  , 


k-l 


(2) 


y^-\/a7  hji)  Pkl  +  y~l  a/o!  ~  1) ^  +  1)  +  ”*:(*)  ! 


/=! 


l=k+l 


1=1 


where  {  nk{i)  }j^— i  is  n  sequence  of  colored  Gaussian  noise  due  to  the  multi-shot  matched  filtering;  and 

pT 

Pki=  Sk{t)si{t-Ti  +  Tk)dt; 

Jo 

^  \  Jo  ~ 

l<k 


(-)  ^ 
Pik  = 


0 

/  Sk(t)  si{t  -  T  -  Ti  +  Tk)  dt ,  k>l 

Jo 


0. 


k  <  I 


After  stacking  the  outputs  of  a  bank  of  K  multi-shot  matched  filters  into  a  /t"  X  1  vector  x(j)  = 

a;i(i)  X2{i)  •••  a;K(i) 

matched  filter  outputs. 


,  we  obtain  the  following  matrix  decomposition  notation  for  the  ith  bit  multi-shot 

x(i)  =  P  A  b(z)  -h  Pc;  A  b(i  -  1)  -1-  Pi  A  b(i  +  1)  +  n(i)  , 
z{i  -  1) 

z(i)  1  +  n(i)  , 
z{i  -1-  1) 


(3) 


where  z(z)  =  Ab(i)  ;  n(i)  ~  A/'(0,  er^  P)  is  a  colored  Gaussian  noise  vector;  P  is  a  symmetric  and  positive  definite 
matrix  defined  by  F{k,  1)  =  pki  ;  Pu  is  an  upper  triangle  matrix  defined  by  Fu{k,  0  =  p\k^ and  Pi  is  a  lower 
triangle  matrix  defined  by  Pi(Ar,  1)  =  p\k^ ■  We  can  even  show  that  Pc;  =  P^- 

Since  the  matrix  A  in  (3),  defined  as  A  =  diag  {  ■■■,  y/a^  } ,  is  a  positive  definite  diagonal  matrix, 

as  long  as  the  sequence  of  vector  {z(z)}  can  be  estimated  properly,  the  multiuser  information  bit  sequence  {b(i)} 
can  then  be  obtained  from  the  following  decision  rule, 


b(i)  =  sign{z(i)}  , 

where  we  assumed  the  multiuser  information  sequence  {b(i)}  is  generated  by  BPSK  source. 

Note  from  (3)  that  for  the  special  case  of  synchronous  channel,  where  the  delays  of  all  K  users  are  equal,  the 
matrices  Pc;  and  P^  vanish.  Hence,  for  the  case  of  synchronous  channel,  only  the  multiuser  interferences  appear 
in  the  signal  term  of  data  x(i).  But  for  the  more  gerenal  case  of  asynchronous  channel,  besides  the  multiuser 
interferences,  intersymbol  interferences  (ISI)  also  exist  due  to  the  additional  terms  in  (3).  Our  problem  is  to 
design  a  near-far  resistant  multiuser  detector  of  reduced  computational  complexity  (linear  in  the  number  of  users 
K)  so  that  we  can  detect  and  separate  multiuser  information  bit  sequence  {b(z)}  from  the  data  sequence  {x(z)}. 
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In  the  next  section,  we  further  analyze  the  structure  inherent  in  the  properly  arranged  multi-shot  matched 
filtering  data,  the  statistics  of  the  noise  component,  and  the  bias  introduced  by  matrix  approximation.  We  then 
propose  a  computationally  simple  approach  to  design  a  near-far  resistant  multiuser  detector  for  asychronous 
CDMA  system. 


Ill  Proposed  Multiuser  Separation  Approaches 


Based  on  the  multi-shot  matched  filtering  output  data  model  in  (3),  we  further  arrange  the  available  N-hit  data 
sequence  {  x(i)  into  a,  N  K  x  1  vector  as  follows, 


x(i) 

x(i  +  1) 
x(z  -I-  2) 

c(i  -f-  AT  -  1) 


=  0  0 


Pi  0  0 

P  Pi  0 


0  0  0 


Pf/  p  Pi 


NKx(N+2)K 


-  1) 

z{i) 

z{i  +  1) 
z{i  +  2) 

z{i  +  n) 


n(z) 
n(z-f  1) 
n(z  +  2) 

n(z  -|-  JV  —  1) 


Pl  0 
P  Pl 
Pc;  P 


z(z) 

z{i  -b  1) 
z{i  +  2) 


Pu  z{i  -  1) 


n(0 

n(z  +  1) 
n(z  +  2) 


Pu  P 


z{i  -b  iV  —  1) 


Pl  z{i  -b  n) 


n(z  -b  AT  —  1) 
noise  N(:,i) 


Note  that  the  first  term  in  (5)  is  an  approximation  of  the  first  term  in  (4).  We  have  proven  that  the  matrix 
V  in  (5)  is  an  (NK  x  NK)  symmetric,  positive  definite,  block-tridiagona,!,  and  full-rank  matrix.  If  we  neglect 
the  bias  term  in  (5),  the  approximation,  a  simple  truncation,  makes  estimation  easy.  Even  further,  we  have 
shown  that  the  noise  vector  in  (5)  is  colored  Gaussian  with  a  probability  density  function  (PDF)  of  Af{0,  TZ), 
and  TZ  —  a'^V.  Hence,  the  Unconstrained  Maximum  Likelihood  Estimate  (MLE)  of  A-bit  multiuser  information 
Z(:,  i)  and  B(:,  i)  are  obtained  from, 


Z(;,  i)  =  (V^  R-i  P)  ^  R-i  X(:,  i)  =  P"!  X(:,  i) , 
B(:,  i)  =  sign|z(:,  z')j  . 


(*■=  1,  2,  3,  ...) 


It  should  be  pointed  out  that  from  (3)  we  can  see  the  ISI  due  to  the  asynchronous  transmission  comes  from  the 
neighboring  two  bits.  That  is,  each  bit  of  multiuser  information  is  present  at  most  in  three  bits  of  multi-shot 
matched  filtering  data.  In  order  to  use  (6)  to  detect  and  separate  multiuser  information  bits,  we  need  to  jointly 
process  M-bits  of  data  x(i),  and  the  minimum  value  of  N  should  be  equal  to  three.  We  further  notice  that  there 
exists  a  nice  structure  in  the  properly  arranged  iV-bit  multiuser  information  |z(:,  z)|  ,  which  can  be  used  to 
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further  improve  the  overall  detection  performance.  We  can  arrange  the  estimates  in  (6)  into  a  square  matrix,  i.  e. 

■  z(l)  z(2)  z(iv) 

..  z(2)  z(3)  •  •  •  z{n  + 1) 

Z=  Z(:,  1)  Z(:,  2)  Z(:,  tv)  J  =  .  .  . 

_  z(iv)  z(jv  +  l)  z{2n) 

The  structure  of  Z  in  (7)  suggests  that  the  improved  estimate  of  the  multiuser  information  bit  b(i)  can  be  obtained 
via  an  anti-diagonal  averaging  processing.  Or, 

b(i)  =  sign  I  ^  ^  Z  (  (tv  —  i)  Jf  +  1  ;  (tv  -  j  +  1)  tv,  o|^  • 

The  simulation  results  in  figure  1  of  next  section  show  the  improvement  in  the  detection  performance  (m  terms 
of  error  probability)  obtained  via  anti-diagonal  averaging. 

Note  that  the  above  results  in  (6),  (7),  and  (8)  are  obtained  by  neglecting  the  bias  term  (caused  by  truncation 
approximation)  in  (5).  Certainly,  the  effect  of  this  truncation  will  be  small  when  we  jointly  process  a  large  number 
of  multiple  bits  of  multi-shot  data  (using  a  large  value  of  N).  Our  simulation  results  in  figure  2  indeed  show 
the  performance  improvement  by  incresing  N .  But  incresing  the  data  dimension  N  will  unavoidably  increase  the 
computations  involved.  From  figure  2  we  also  observe  that  there  exits  a  limit,  which  is  near-far  resistant,  on  the 
performance  even  though  we  can  keep  increasing  the  value  of  N .  And  this  limit  is  due  to  the  existence  of  noise 
term  in  (5).  Our  observation  is  that  with  the  help  of  the  inherent  structure  in  matrix  V,  we  can  handle  the  bias 
term  properly,  and  still  achieve  the  near-far  resistant  performance  limit  without  increasing  the  value  of  N  and 
computations  involved.  Specifically,  matrix  "P  in  (5)  has  a  block-tridiagonal  structure.  Since  P  is  a  symmetric 
positive  definite  matrix,  is  a  lower  triangle  matrix,  and  Pu  is  an  upper  triangle  matrix  with  Py  =  P^^.  We 
can  show  that  the  above  matrix  P  is  a  symmetric  positive  definite  matrix  of  full  rank.  Of  most  importance  is  the 
fact  that  this  matrix  P  is  also  a  diagonal  dominant  matrix.  Therefore,  its  inverse  P~^  is  also  a  diagonal  dominant 
matrix  (see  figure  4).  We  also  notice  the  fact  that  the  bias  vector  in  (5)  is  almost  a  zero  vector,  except  for  the 
first  sub-block  and  the  last  sub-block.  With  the  inverse  matrix  filtering  of  (6),  the  estimate  of  TV-bit  multiuser 
information  has  the  form  of, 

z{i)  P(7z(i  — 1) 

z(i  1)  o' 

Z(:,  f)  =P“^X(:,  i)  =  i  -1-P“^  -fP-’^  N(:,  i)  .  (9) 

. 

z(i  4-  JV  -  2)  0  noise 

z(i-|-JV  — 1)  _  _  Pjrz(i-I-Ar) 

' - ^  ' - - - ^ 

true  value  bias 

Due  to  the  diagonal  dominant  feature  of  matrix  P~\  the  estimates  of  the  first  bit  multiuser  information,  z(f), 
and  the  last  bit  multiuser  information,  z{i  -f  JV  —  1),  will  be  mostly  affected  by  the  truncation  bias.  While  the 
estimate  of  multiuser  information  bit  in  the  middle  sub-block  of  the  TV-bit  vector  Z(;,  i)  will  be  least  affected  by 
the  truncation  bias.  Therefore,  we  can  simply  use  the  minimum  value  of  TV  =  3,  and  only  estimate  the  middle 
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bit  of  multiuser  information  with  reduced  computation  as  follows.  We  denote  the  inverse  of  T’  as 


'  p 

Pi 

0 

-1 

‘All 

Ai  2 

Ai3 

Pu 

P 

Pl 

_ 

A21 

A22 

A23 

0 

Pu 

P 

.  A31 

A3  2 

A33 

(10) 


Since  we  need  only  to  estimate  the  middle  bit  of  multiuser  information  z(z  +  1)  by  joint  processing  of  3-bit 
multi-shot  matched  filtered  data  vector  x(i),  x(i  -f  1),  x(i  -|-  2),  i.  e. 


In  general,  the  K  X  NK  matrix  W  can  be  calculated  in  advance  by  iteratively  using  the  results  of  block  matrix 
(of  less  dimension)  inversion  and  the  nice  structure  of  matrix  V  in  (5).  t'or  the  case  of  iV  =  3,  we  calculated 
K  X  3Jf  matrix  W  as  follows, 


W  = -(P -PiP  iP[/ -  [  Pt;P-i  -I  p^p-i  j  .  (12) 

In  figure  5,  we  demonstrate  the  system  structure  of  the  signal  processing  parts,  base  on  (11)  and  the  matrix  W  in 
(12).  It  can  be  seen  that  when  the  transmission  delay  changes,  instead  of  re-inverting  the  whole  (NKxNK)  matrix 
V ,  we  only  need  to  update  the  sub-matrices  P,  Py,  Pjj  of  less  dimension  with  much  less  computations.  Simulation 
results  in  figure  3  also  show  the  near-far  resistant  performance  achieved  by  the  proposed  computationally  simple 
multiuser  detector  with  debiasing  technique  of  (11). 


IV  Simulation  Results 

Simulation  results  using  the  proposed  multi-shot  approaches  are  shown  in  the  figures  1  through  3.  In  these  figures, 
we  compare  the  detection  performance  of  the  proposed  detector  with  that  of  the  conventional  detector.  We  also 
investigate  the  potential  improvement  in  performance  with  further  anti-diagonal  averaging  processing;  the  effects 
of  increasing  the  matrix  dimention  (value  of  N)',  and  the  effects  of  debiasing  the  approximation  error  to  get  a 
near-far  resistant  dectector.  In  our  simulations,  we  choose  normalized  Gold  code  of  length  seven  as  the  normalized 
signature  waveforms  for  the  case  of  three  users, 

1  1 

Si  =  ;^  [  1  -1  -1  1  1  1  -1  J  ; 

1  ‘  ‘  1 ' 

The  transmission  delays  associated  with  all  three  users  are  chosen  as, 

n  =  0,  r2  =  Tc,  T3  =  2T(, , 
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where  Tc  stands  for  the  chip  interval,  and  the  symbol  interval  T  =  7%  for  our  simulations. 

Note  that  the  SNR  of  each  user  is  defined  as  the  ratio  of  bit  energy  to  noise  density,  or, 

SNR{k)  =  ^,  (i  =  l,  2,  3) 

(T 

where  we  assumed  the  symbol  interval  T  is  normalized. 

In  each  example,  100, 000  independent  trials  are  run  for  obtaining  the  error  probability  of  the  desired  user.  Shown 
in  figure  1  through  figure  3  are  the  error  probability  of  the  desired  user  (with  fixed  SNR)  versus  the  SNRs  of  each 
interferencing  user,  with  various  multiuser  processor.  Also  shown  in  these  figures  are  the  detection  performance 
of  the  conventional  single  user  detector  for  comparision.  The  near-far  resistant  feature  of  the  proposed  multiuser 
with  reduced  computation  detector  can  be  easily  seen  from  figure  3. 

V  Conclusions 

In  this  paper,  we  proposed  a  computationally  simpe  multiuser  detection  and  separation  scheme  for  asynchronous 
CDMA  communication  applications.  The  proposed  detector  achieves  the  near-far  resistant  performance  with  a 
very  simple  debiasing  signal  processing  technique.  The  possible  simplifica,tion  in  system  implementation  is  also 
demonstrated  by  incorporating  the  block-tridiagonal  structure  of  the  matrix  V  into  our  debiasing  technique. 
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Pe(1) 


Figure  1:  Performance  comparison  of  various  detectors  in  asynchronous  multiuser  CDMA  communication  systems. 
The  performance  improvement  obtained  by  anti-diagonal  averaging  processing  is  also  shown  in  this  Sgure. 
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Pe(1)  versus  SNR(i)-SNR1  in  dB 
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Figure  2:  Performance  comparison  of  various  detectors  in  asynchronous  multiuser  CDMA  communication  systems 
The  performance  improvement  obtained  by  increasing  the  data  dimension  N  is  shown  in  this  figure. 


Pe(1)  versus  SNR(i)-SNR1  in  dB 
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Figure  3:  Performance  comparison  of  various  detectors  in  asynchronous  multi-user  CDMA  communication  sys 
terns.  The  performance  improvement  obtained  by  debiasing  processing  is  shown  in  this  figure. 
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(b).  structure  of  matrix 


Figure  4:  Structure  of  matrix  P  and  matrix  P  ^  shown  by  mesh  and  contour  plots.  Parameters  used  are  the 
same  as  that  used  in  figure  1.  Note  that  the  diagonal-dominant  feature  appears  in  both  V  andp-^. 
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